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Foreword  (0) 
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(U)  The  Acoustic  Model  Evaluation  Committee  (AMEC)  haa 
been  chartered  to  serve  as  an  advisory  group  to  the 
Director,  Naval  Oceanography  Division  (OP-952),  on 
matters  dealing  with  model  evaluation.  In  fulfillment 
of  it's  charter,  AMEC  will  produce  a series  of  reports 
detailing  the  results  of  model  evaluations.  Volume  I 
described  the  evaluation  methodology  selected  and  the 
manner  in  which  it  lias  been  Implemented.  Volume  IA 
describes  experimental  propagation  loss  data  sets 
suitable  for  the  evaluation  of  models  in  a range 
dependent  environment.  Volume  II  presented  the  results 
of  evaluating  the  PACT  PL9D  propagation  loss  model. 
This  report.  Volume  III,  presents  the  results  of 
evaluating  the  RAYWODE  X propagation  loos  model. 

. O.  -Qq. 


G.  T.  Phelps,  Captain,  USN 
Gonriding  Officer 
NGSIDA 
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(U)  The  Acoustic  Model  Evaluation  Committee  (AMEC)  has 
applied  the  methodology  described  in  Volume  I of  this 
series  of  reports  to  evaluate  the  RAYMODE  X propaga- 
tion loss  model.  The  accuracy  of  RAYMODE  X has  been 
assessed  by  quantitative  comparisons  with  eight  set  of 
experimental  data  covering  a broad  spectrum  of  en- 
vironmental acoustic  scenarios.  The  physics  of  RAYMODE 
X has  been  examined  by  R.  Deavenport  of  the  Naval 
Underwater  Systems  Center,  New  London  Laboratory. 
RAYMODE  X typically  has  run  times  between  5 and  30 
seconds  on  the  UNIVAC  1108  computer.  The  model  is 
poorly  documented  with  the  exception  of  a well-written 
user's  guyide;  this  extends  to  a severe  lack  of  com- 
ment cards  within  the  computer  code.  The  program  could 
clearly  benefit  from  an  improved  surfrce  duct  module; 
no  other  serious  deficiencies  in  che  physics  of 
RAYMODE  X have  been  noted.  Various  versions  of  RAYMODE 
exist  in  fleet  operations:  These  versions  do  not  con- 
tain identical  physics,  are  written  in  different 
computer  languages,  and  are  fun  on  hardware  utilizing 
different  word  lengths.  Consistency  of  results  from 
these  versions  has  not  been  demonstrated.  It  is 
recommended  that  a program  of  configuration  management 
be  aplied  to  all  RAYMODE  versions.  RAYMODE  X has  many 
useful  features  including  eigenray  information,  inde- 
pendence of  Initial  range  and  range  increment  for 
propagation  loss  calculations,  provision  for  vertical 
beampatterns  and  the  ability  to  input  an  external 
bottom  loss  table.  This  evaluation  was  completed  in 
September  1980. 
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Preface  (U) 


(U)  This  report  was  prepared  under  the  joint  sponsor-  * 

ship  of  the  Naval  Sea  Systems  Command,  Program  «. 

Manager,  P.  R.  Tiedeman  (SEA  63D3),  PE  63708N;  the 
Surveillance  Environmental  Acous  ic  Support  Project, 

Program  Manager,  Dr.  Robert  A.  Gardner  (NORM  Code 
520),  PE  63795N;  the  Tactical  ASW  Environmental 
Acoustic  Support  Project,  Program  Manager,  E.  D. 

Chaika  (NORM  Code  530),  PE  63795N;  via  the  auspices  \ 

of  OP-952D  (Capt.  J.  Harlett) . : 
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The  Acoustic  Model  Evaluation  Committee  (AMEC)  Reports 

Volume  III,  Evskiation  of  the  RAYMOBE  X Promotion  Loss  Model  (U) 


1.0  (U)  Introduction  submarine  fleet  and  is  the  transmission 


(U)  This  volume  is  third  in  a series  of 
Acoustic  Model  Evaluation  Committee 
(AMEC)  reports.  Volume  I deals  in  detail 
with  the  model  evaluation  methodology 
and  its  implementation  in  fulfillment  of 
AMEC's  charter.  This  volume  details  the 
application  of  that  methodology  for  the 
evaluation  of  the  RAYMODE  X model  as  run 
on  the  UNIVAC  1108  computer  at  the  Naval 
Underwater  Systems  Center,  New  London 
Laboratory.  The  RAYMODE  X evaluation  was 
completed  on  30  September  1980.  No  modi- 
fications of  RAYMODE  X were  required  in 
order  to  perform  the  accuracy  assessment 
portion  of  the  evaluation;  of  Importance 
in  this  regard,  the  RAYMOCE  X model  has 
the  capability  of  accepting  an  external 
bottom  los  s table  input  and  provision 
for  writing  propagation  loss  versus 
range  results  to  an  external  file. 

^ IT  ^ Tho  mml  nof-T  Art  mnt-KA/Irtl  nnu  -la 
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described  in  Section  1.1  of  this  volume 
and  in  greater  detail  in  Volume  I of 
this  series.  The  primary  areas  in  which 
we  seek  to  provide  model  evaluation 
information  are  (1)  model  description, 
(2)  physics  and  mathematics,  (3)  run 
time,  (4)  core  storage,  (5)  complexity 
of  program  execution,  (6)  ease  of  ef- 
fecting program  alterations,  (7)  ease  of 
implementation  (on  a different  compu- 
ter), (8)  cognizant  individual^)  or 
organizational  el  ment(s),  (9)  byprod- 
ucts, (10)  special  features,  (11)  ref- 
erences, and  (12)  accuracy  assessment. 

(U)  RAYMODE  X is  a computer  program  for 
the  prediction  of  transmission  loss  ver- 
sus range  in  an  environment  which  can  be 
characterized  by  a flat  bottom  and  a 
single  sound  speed  profile.  This  model 
is  in  extensive  fleet  usage,  supporting 
the  Optimum  Mode  Selection  systems  for 
TRIDENT,  the  BQQ-5,  and  the  BQQ-6  in  the 


loss  module  in  the  Sonar  In  Situ  Mode 
Assessment  System  (SIMAS)  used  in  the 
surface  fleet.  RAYMODE  X has  been  dis- 
tributed to  other  naval  activities  and 
Navy  contractors. 

(U)  Four  classes  of  models  are  described 
by  Hersey  (1977):  Research  Model,  Candi- 
date Model,  Navy  Evaluated  Model,  and 
Navy  Operational  Model.  With  the  publi- 
cation of  this  report,  RAYMODE  X has 
fulfilled  the  requirements  for  status  as 
a Navy  Evaluated  Model.  The  RAYMODE  mod- 
el (in  a variety  of  versions)  has  been  a 
Navy  Operational  Model  for  many  years. 

(U)  As  we  shall  see  below,  RAYMODE  X 
typically  has  run  times  ranging  from  5 
to  30  seconds  on  the  UNIVAC  1108  with 
the  EXEC  VIII  operating  system.  The  core 
required  by  RAYMODE  X as  dimensioned  for 
400  points  and  50  modes  and  ray  paths  is 

17110  rl  n a T mnl  r.*A  f/l  n f nv  al  itnlitn  a f a1  Af 
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routines).  These  run  times  and  core 
requirements  allow  RAYMODE  to  be  imple- 
mented in  desktop  calculators  (HP  9845 
and  Tektronix  4051)  and  to  meet  opera- 
tional requirements  for  run  time.  As  of 
the  termination  of  this  evaluation  on  30 
September  1980,  RAYMODE  documentation, 
both  internal  and  external  to  the  com- 
puter code,  was  lacking  with  the  excep- 
tion of  a user's  guide  (Yarger,  1976) 
and  a technical  memorandum  by  Leibiger 
(1971)  that  ptesented  the  essence  of  the 
RAYMODE  method  but  did  not  tie  in  to  the 
computer  code.  The  theory  of  RAYMODE  is 
further  described  by  DiNapoli  and 
Deavenport  (1980). 

(U)  The  version  of  RAYMODE  herein  evalu- 
ated is  RAYMODE  X as  resident  and  run  on 
the  UNIVAC  1108  computer.  Other  RAYMODE 
versions  were  not  tested  but  are  descri- 
bed in  section  10.  References  to  these 
versions  are  provided  as  available. 
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1.1  (II)  Thi  AMEC  «iM% 

(U)  Volume  I of  this  series  of  reports 
presents  the  AMEC  model  evaluation 
methodology  in  detail.  The  following 
list  Is  a synopsis  of  this  methodology; 
the  Items  are  taken  from  an  information 
request  form  sent  to  those  persons  re- 
sponsible (usually  the  developer)  for  a 
model  which  Is  to  undergo  evaluation 
and,  taken  together  with  the  physics 
review  and  accuracy  assessment,  consti- 
tute the  evaluation. 

(U)  Range  Independent  Propagation  Loss 
Information  requested  for  AMEC: 

1.  (U)  Model  Description 

• Purpose(s)  of  the  model. 

• List  of  Input  variables  and  their 
units  (Inputs  obtained  from  associated 
data  bases,  internal  routines,  functions 
or  tables  should  be  so  identified). 

• List  of  output  options.  Examples  of 
tabular  and  graphical  results. 

• A list  of  systems  (e.g.,  sonar  predic- 
tion, engagement  model,  etc.)  supported 
by  the  model,  including  the  role  of  the 
model  in  the  system  and  the  stated  pur- 
pose of  the  systems. 

• Limitations  designed  into  the  model, 
through  inherent  limits  of  the  physics, 
mathematics,  environm  ntal  description, 
computer  implementation,  etc.  These  lim- 
itations, taken  together,  define  the 
model' 8 domain  of  applicability  and  in- 
clude frequency,  bandwidth,  range,  etc. 
Also  included  are  limitations  involving 
choice  of  computer,  graphics,  and  telem- 
etry links.  Please  outline  the  extent  to 
which  the  limitations  result  from  design 
decisions  based  upon  the  basic  purpose 
of  the  model  development  effort  and/or 
tradeoffs  required  by  time  (run  time  or 
product  delivery),  cost,  and  computer 
assets . 

• A list  of  extant  model  versions.  Note 
differences  between  versions  including 
computer,  changes  in  inputs  and  outputs, 
assignments  of  default  values,  graphics, 
program  language,  use  of  overlays,  etc. 


2.  (U)  Run  Time 

• Provide  run  time  as  function  of  com- 
puter, number  of  points  and  inpu t/output 
selections,  and  model  version.  Divide 
run  time  into  computation  time  and  time 
required  for  printing  and  plotting. 
Describe  tradeoffs  ;tween  accuracy  and 
run  time  as  affected  by  input  options. 

3.  (U)  Core  Storage 

• Provide  information  on  core  storage 
requirements  on  a version  basis.  Identi- 
fy techniques  used  to  reduce  core 
requirements  including  use  of  overlays, 
memory  mapping,  disk  memory  swap  and  the 
use  of  techniques  such  as  interpolation 
in  place  of  calculation. 

4.  (U)  Complexity  of  Execution 

• Provide  a program  listing. 

• Define  all  input  and  output  parameters 
under  user  control. 

• What  default  values  or  conditions  are 
assigned  with  the  program? 

• Identify  restrictions  on  parameter 
values . 

• Identify  unusual  parameters  and  pro- 
vide guidance  for  their  selection. 

• Does  a user's  guide  exist?  If  so, 
please  give  reference. 

5.  (U)  Ease  of  Effecting  Program  Alterations 

• Supply  a program  flow  chart. 

• A list  of  program  variables  and  their 
definition. 

• Extent  to  which  a model  is  tied  into  a 
specific  computer  executive  system  or 
special  equipments  or  programs,  library 
routines,  etc. 

6.  (U)  Ease  of  Implementation  on  a Different 
Computer 

• List  of  computers  (and  executive  sys- 
tems) on  which  model  is  presently  run- 
ning . 

• List  of  military  and  civilian  activi- 
ties using  the  model. 
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• Computer  language(s)  used  by  the  model 
(all  versions). 

• Special  co  es  (e.g.,  plotting  rou- 
tines, library  functions). 

• To  what  extent  is  the  program  depend- 
ent on  a given  computer  executive  sys- 
tem? 

• Identify  test  cases  to  assure  proper 
running  on  a new  computer  (including 
scenarios  treated);  are  all  subroutines 
and  lines  of  code  exercised. 

• A list  of  all  errors  returned  and  the 
situations  which  cause  them. 

7,  (U)  Cognizant  Individual(s)  and/or 
Organizational  Elements,  Names  and  Addresses 
of  Those  Reponsible  for 

• Theory  upon  which  model  is  based. 

• Model  development. 

• Computer  Implementation. 

• Model  maintenance  and  conf iguration 
management . 

8.  (U)  References 

(U)  A list  of  references,  including 
those  which  discuss  theories  upon  which 
model  is  based,  and  numerical  methods 
employed.  References  worthy  of  special 
mention  are  (a)  a user’s  guide;  and  (b) 
a response  to  SECNAVINST  3560.1,  Tac- 
tical Digital  Systems  Documentation 
Standards  of  8 August  1974,  or  to  DOD 
Standard  7935. 1-S,  Automated  Data  Sys- 
tems Documentation  Standards  of  13  Sep- 
tember 1977,  or  to  other  Navy  or  DOD 
documentation  requirements. 


9.  (U)  By-Products 

• A list  of  output  by-products  (e.g., 
eigenray  information,  arrival  angle  vs. 
range,  ray  diagram). 

• A list  of  by-products  not  available 
externally  but  which  are  internally  cal- 
culated . 

10.  (U)  Special  Features 

• A list  of  special  features  (e.g.,  pro- 
vision for  beampatterns , multi-frequency 
results  through  interpolation,  etc.). 


(U)  The  review  of  the  physics  and  math- 
ematics and  computer  implementation  of  a 
given  model  i9  undertaken  by  an  Inde- 
pendent expert  in  the  appropriate  field 
of  modeling.  In  particular,  the  physics 
and  mathematics  are  examined  to  define 
the  model's  domain  of  applicability 
throvtgh  assumptions,  approximations  and 
the  assignment  of  "nominal  values"  to 
various  parameters. 

(U)  The  reporting  of  the  model's  physics 
includes  the  basic  foundations  and  ap- 
proach and  any  unusual  techniques  and, 
especially,  any  extensions  to  theory 
and/or  unique  capabilities  otherwise 
unavailable.  Examination  of  the  model's 
physics  and  mathematics  is  to  include 
consideration  of  environmental  inputs, 
including  theories  and  the  appropriate- 
ness of  data  base  selection.  The  com- 
puter implementation  is  examined  to 
assure  that  the  calculations  required  by 
tht  theory  are  correctly  performed.  The 
efficiency  or  other  aspects  of  the  pro- 
gram code  are  not  addressed  here. 

(U)  Two  accuracy  assessment  procedures 
are  employed  in  AMEC  evaluation.  Both 
yield  quantitative  results  and  involve 
comparison  of  model  outputs  with  experi- 
mental data  or  the  output  of  a reference 
model.  The  steps  of  the  first  procedure, 
called  the  Difference  technique,  follow: 
(1)  Smooth  the  reference  data  set  (only 
if  CW  or  exhibiting  large  fluctuations) 
and  the  output  of  the  model  (only  if 
coherent  phase  addition  was  used)  by 
applying  a 2 kra  moving  window.  (2)  Sub- 
tract the  model  output  from  the  refer- 
ence data  set  (after  appropriate  smooth- 
ing). (3)  If  possible,  divide  the  dif- 
ference curve  into  range  intervals  cor- 
responding to  direct  path,  bottom  inter- 
action and  convergence  zone  modes.  If 
not  possible,  either  (a)  do  not  subdi- 
vide into  range  intervals;  (b)  use 
quasiarbitrary  intervals,  which  may  be 
tactically  useful;  or  (c)  subdivide  on 
the  basis  of  any  features  evident  in  the 
measured  data.  (4)  In  each  range  inter- 
val calculate  the  mean  ^ and  the  stand- 
ard deviation  a of  the  differences.  (5) 
Analyze  results,  attempting  to  identify 
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causes  of  discrepancies.  The  above  steps 
are  supported  by  figures  as  follows: 
measured  data,  smoothed  measured  data, 
model  output,  smoothed  model  output,  and 
difference  between  smoothed  curves. 
These  curves  are  drawn  to  the  same  scale 
and  may  be  overlaid  on  a light  table, 
facilitating  eyeball  comparison  and 
diagnosis.  As  useful  as  this  technique 
is  in  identifying  significant  differ- 
ences and  facilitating  diagnosis,  it  has 
a number  of  shortcomings:  (a)  misleading 
in  convergence  zones  where  range  errors 
are  as  significant  as  errors  in  level; 
(b)  it  is  conceivable  that  large  errors 
occur  at  dB  levels  of  no  consequence  for 
operational  systems;  and  (c)  the  differ- 
ence approach  leads  to  answers  which  are 
not  particularly  useful  to  fleet  pur- 
poses, especially  in  the  context  of  spe- 
cific sonar  systems.  These  shortcomings 
are  eliminated  in  the  second  accuracy 
assessment  technique,  called  the  FOM 
(Figure  of  Merit)  technique.  In  this 
technique  the  data  is  once  again  smooth- 
ed as  in  step  (1)  above.  FOM  are  then 
selected  in  5 dB  steps.  For  each  FOM, 
detection  range  Information  is  tabu- 
lated: range  of  continuous  coverage, 
ranges  of  convergence  zone  starts  and 
ends,  and  in  range  intervals  over  which 
detection  coverage  is  zonal  in  nature — 
the  percentage  of  the  interval  over 
which  detection  can  be  made.  This  FOM 
vs.  detection  range  analysis  is  perform- 
ed for  model  and  reference  data  set,  the 
results  compared  and  reasons  sought  for 
significant  disparities. 

(U)  Taken  together,  the  two  accuracy  as- 
sessment techniques,  the  Difference  and 
FOM  techniques,  lead  to  results  useful 
to  scientific  analysis  and  for  system 
performance  estimation. 

2.0  (U)  RAYMODE  X Description 

(U)  RAYMODE  X is  the  most  recent  in  an 
evolutionary  chain  of  RAYMODE  models  and 
is  the  basic  version.  The  terra,  RAYMODE, 
indicates  a method  of  calculating  propa- 
gation loss  that  utilizes  both  ray  and 
normal  mode  theories.  In  addition  to  the 
calculation  of  propagation  loss  versus 


range,  RAYMODE  also  calculates  arrival 
angle  versus  range  and  travel  time 
versus  range  for  the  various  rays  and 
groups  by  index:  surface  duct  (J"l), 
convergence  zone  (J“2),  and  bottom 
bounce  (J"3)  paths. 

(U)  The  model  is  dimensioned  to  give  a 
maximum  of  400  range  points  at  which 
propagation  loss  is  calculated  (using 
either  coherent  or  incoherent  phase  ad- 
dition). These  range  points  are  equally 
spaced  between  user-specified  minimum 
and  maximum  ranges.  Beam  patterns  may  be 
input  for  both  source  and  receiver.  Al- 
though RAYMODE  has  an  internal  subrou- 
tine from  which  bottom  loss  versus  graz- 
ing angle  is  specified  given  a bottom 
type  (1  to  9),  a bottom  loss  versus 
grazing  angle  table  may  be  input  from  an 
external  source.  The  bottom  loss  curves 
were  obtained  from  Marine  Geophysical 
Survey)  (MGS)  data  (Podezwa,  1975).  RAY- 
MODE contains  a subroutine  for  the  cal- 
culation of  surface  loss  (Beckman  and 
Spizzichino,  1963).  As  run  at  NUSC,  RAY- 
MODE contains  graphics  routines  utiliz- 
ing Integrated  Graphics  System  (IGS) 
software  and  the  Information  Interna- 
tional FR80  hardware. 

(U)  A list  of  physical  variables  used  in 
the  RAYMODE  program  is  given  in  Table 
2-1.  The  Fortran  label  (or  variable  is 
given  in  the  left-hand  column  and  the 
definition  in  the  right-hand  column. 

(U)  In  addition  to  a sound  speed  versus 
depth  or  temperature  versus  depth  plus  a 
constant  salinity  value,  the  program  can 
access  historical  sound  speed  data 
fields.  Adjacent  equal  sound  speeds  or 
temperatures  are  modified  to  avoid  a 
zero  gradient  condition. 

(U)  Systems  supported  by  the  versions  of 
the  RAYMODE  model  (not  RAYMODE  X as  it 
exits  on  the  UNIVAC  1108)  are: 

• Optimum  Mode  Selection  (OMS)  for  TRI- 
DENT 

• Improved  Sonar  Processing  Equipment 
( ISPE ) — sonar  suite  replacement  on  SSBNs 
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ACOS: 
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ALOG-, 

AL0G10 : 

</ 

AMAX1: 

k 

AMIN1: 

wm 

AMPTUD: 

85 

B 

ANGLE: 

v< 

a 

ANGLN: 

AEGLX: 

ts 

ANGLO: 

9 

AR: 

BEAM: 

K 

BIGM: 

R 

BL: 

Oft 

BPI: 

1 

C: 

CA: 

r<i 

V; 

te. 

CAD EL: 

t--' 

CAMAX: 

u 

CAMAX -T; 

K' 

CAMIN: 

CAMINJ : 

1 

CAMU: 

CAN: 

Table  2-1.  (U)  Physical  variables  used  in  the  RAYMODE  program  ;'| 

Absolute  value  function  in  Univac  Library;  takes  real  argument,  gives  [I 

real  result . t 

s 

Arccos  function  in  Univac  Library,  -l._<  argument  <1. , result  in  ft 

radians . 1 3 

Natural,  log  function  in  Univac  Library.  ft 

Log--  in  Univac  Library.  jP 

ft 

Univac  function  to  select  maximum  of  2 real  numbers.  ^ 

Univac  function  to  select  minimum  of  2 real  numbers.  ^ 

Subroutine.  P! 

:■] 

+ Maximum  sonar  angle  in  degrees  or  velocity  terms  (input — see  g 

Yarger  (19"6)).  !;3 

+ Minimum  sonar  angle  in  degrees.  j| 

+ Maximum  sonar  angle  in  degrees.  ft 

+ Minimum  sonar  angle  in  degrees  or  velocity  terms  (input — see  ft 

Yarger  (1976)).  ft 

Internally  computed  range  value  in  yards.  w 

ft 

Subroutine.  ft 

Largest  mode  trapped.  ft 

Array  of  bottom  loss  values  in  dB  (may  be  input —see  Yarger  (1976)).  * 

Large  integral  multiple  of  2tt  for  use  in  SQUD  subroutine.  *1 

ft 

fpj 

Array  of  profile  velocities  (input — see  Yarger  (1976)).  jjj 

Intermediate  wave  number.  1 

ft 

ft 

Difference  between  wave  numbers.  ft 

ft 

Maximum  wave  number  trapped  by  an  angular  interval.  n] 

t 

Array  to  save  CAMAX  values  to  provide  plot  routines.  ^ 

Minimum  wave  number  trapped  by  an  angular  interval.  y 

V 

Array  to  save  CAMIN  values  to  provide  plot  routines.  ft 

i 

Wave  n mber  associated  with  receiver  depth.  C 

hJ 

Wave  number  associated  with  bottom  depth. 
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CAONE: 


CMAX: 


CMIN: 


DELIH: 


DELRE: 


DELTAR: 


DJCUT: 


Wave  number  associated  with  source  depth. 

Wave  number  associated  with  surface  depth. 

Wave  number  associated  with  maximum  sonar  angle. 

Wave  number  values  defined  at  points  along  a path. 

Wave  number  associated  with  velocity  CQ. 

Maximum  velocity  trapped  by  an  angular  interval. 

Minimum  velocity  trapped  by  an  angular  interval. 

Slightly  adjusted  velocity  to  search  for  phase  changes. 

Univac  Library  function  for  trigonometric  cosine;  argument  in 
radians . 

Array  of  adjusted  profile  velocities  in  yards  used  internally. 

Maximum  velocity  between  source  and  receiver,  slightly  shifted. 

Term  containing  surface,  bottom,  and  beam  losses  for  bottom  bounce. 

Range  derivative  array  associated  with  cycle  range  at  each  point 
along  a ray. 

Range  derivative  array  associated  with  range  adjustment  to  cycle 
range  for  receiver  depth. 

Constant  number  of  degrees  per  radian. 

Imaginary  part  of  complex  propagation  loss  contribution. 

Real  part  of  complex  propagation  loss  contribution. 

Range  increment  (input-see  Yarger  (1976)). 

Array  of  terms  containing  amplitude,  surface  and  bottom  losses  for 
bottom  bounce. 

Cutoff  to  determine  if  losses  are  so  high  as  to  be  impractical  to 
compute. 

Array  of  losses  in  dB  for  source  deviation  loss  table  (input — see 
Yarger  (1976)). 

Computed  beam  loss. 

Array  of  losses  in  dB  for  receiver  deviation  loss  table  (input — see 
Yarger  (1976)). 
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DS: 


DTAPLT: 

EMDEL: 

ENL14: 

EXITG: 

EXP: 

F: 

FIRST: 

FLOAT: 

G: 

GETPRO: 

HALFPI: 

HC: 

HD: 


HDP: 

HDR: 

HS: 


HSP: 

IDL: 

IFIX: 


IFLAG: 


IJDL: 


INPUTS: 

INT: 


Range  derivative  array  associated  with  range  adjustment  to  cycle 
range  for  source  depth. 


Subroutine. 

Difference  between  mode  numbers. 

Factor  determining  spacing  between  wave  number  points  along  a path. 
IGS  subroutine  to  terminate  plotter. 

Univac  Library  for  to  power. 

Frequency  in  Hz  (input — see  Yarger  (1976)). 

Subroutine  to  Initialize  time  clock  to  0. 

Univac  Library  function  to  convert  integer  to  real  variable. 
Gradient . 

Subroutine. 

tt/2. 

Phase  array  associated  with  cycle  range  points  along  a path. 


Phase  array  for  receiver  depth  adjustment  associated  with  cycle 
range  points  along  a path. 


Point  linearly  interpolated  from  HD  array. 

Header  array  in  format  12A6  (input — see  Yarger  (1976)). 


Phase  array  for  source  depth  adjustment  associated  with  cycle  range 
points  along  a path. 


Point  linearly  interpolated  from  HS  array. 

Number  of  points  in  source  beam  pattern  (input — see  larger  (1976)). 


Univac  Library  function  to  convert  real  numbers  to  integer,  i.e., 
truncated. 


Flag  set  by  subroutine  LINTRP  Indicated  point  to  be  interpolated 
already  existed. 


Indicator  if  >o  that  either  source  or  receiver  (or  both)  beam 
patterns  exist. 


Name  of  NAMELIST  inputs  set. 

Univac  Library  function  to  convert  real  argument  to  integer. 
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IOC SAN: 
IPRINT: 
IPROFL: 

IQ: 

IQMAX: 

IQMIN: 

IQO: 

IR: 

ISEASN: 

ITAB: 


J: 

JDL: 

K: 

KEEPMU: 

KEEPNU: 

L: 

LAMDA: 

LAMDAB: 

LAMDAJ: 

LAMMIN: 

LEROY: 

LINTRP: 

H: 

MAX: 

MAXJ: 


Ocean  code  (Input — see  Yarger  (1976)). 

Print  option  (input — see  Yarger  (1976)). 

Profile  code  (input — see  Yargor  (1976)). 

Cycle  number  index. 

Maximum  cycle  number. 

Minimum  cycle  number. 

Cycle  number  internal  value. 

Range  index. 

Season  code  (input — see  Yarger  (1976)). 

Number  of  points  in  bottom  loss  table  (roavbe  input — see  Yarger 
(1976)). 

Index  of  angular  interval  trapped  by  profile. 

Number  of  points  in  receiver  beam  pat' er  w 

Generally  used  to  index  points  along  a iath  or  wave  numbers. 

Save  receiver  index  MU  as  input. 

Save  source  index  NU  as  input. 

Path  index. 

Number  of  cycles  for  non-bottom  bounce  (input — see  Yarger  (1976)). 

Number  of  cycles  for  bottom  bounce  (input — see  Yarger  (1976)). 

Array  to  save  number  of  cycles  for  each  angular  interval  to  pass  to 
ray  and  travel  time  plot  routines. 

Minimum  cycle  (input — see  Yarger  (1976)). 

Subroutine. 

Subroutine. 

Mode  index. 

Univac  Library  function  to  pick  maximum  of  2 integers. 

Parameter  to  show  maximum  number  of  allowable  angular  intervals  for 
dimensions. 
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am 

MAXK: 

I 

MAXM: 

1 

MAXMOD: 

MAXN: 

MAXR: 

i 

MBIGST: 

it« 

METRIC: 

g 

MGSBL : 

MG SOP: 

1 

MINMOD: 

Li] 

MJ: 

S3 

■i 

MM1: 

I 

MODESG: 

MS: 

a 

MSMLST: 

| 

MU: 

MO: 

i 

N: 

i v 

NEGB: 

ft 

NKEXP: 

A 

r>l 

u 

NL: 

.•  ■ 

NLP: 

if? 

< 
!*'  - 

NLPJ: 

g 

NMODE: 

C;.j 

NP: 

Parameter  to  show  maximum  number  of  points  along  a path  for 
dimensions. 

Parameter  to  show  maximum  number  of  modes  trapped  for  dimensions. 
Maximum  mode  number  (input — see  Yarger  (1976)). 

Parameter  to  show  maximum  number  of  allowable  points  in  profile  and 
input  tables. 

Parameter  to  show  maximum  number  of  ranges  allowed  for  dimensions. 
Biggest  mode. 

Metric  input/output  option  (input — see  Yarger  (1976)). 

Subroutine. 

MGS  bottom  loss  province  (input — see  Yarger  (1976)). 

Minimum  mode  number  (input— see  Yarger  (1976)). 

Number  of  modes  between  biggest  and  smallest. 

Number  of  bottom  bounce  J index. 

IGS  subroutine  to  Initialize  plotter. 

Number  of  points  in  surface  loss  table. 

Smallest  mode  number. 

Receiver  index  on  profile  (may  be  input — see  Yarger  (1976)). 

Mode  cutoff  (input — see  Yarger  (1976)). 

Number  of  points  in  input  profile  (input — see  Yarger  (1976)). 

Indicator  to  sense  numbers  beyond  certain  limits  for  mode  summation. 

Exponent  to  determine  spacing  of  wave  numbers  along  a path  (non- 
bottom) . 

Number  of  points  along  a path. 

Adjusted  number  of  points  along  a path. 

Array  to  save  NLP  per  angular  interval  to  pass  to  ray  and  travel 
time  routines. 

Number  of  modes. 

Adjusted  number  of  points  in  velocity  profile. 
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NR: 

NU: 

NXPO: 

OMEGA: 

CNE3RD : 

PHI: 

PHI1: 

PHI2: 

PI: 

PL: 

PLIM: 

PLOTCZ: 

PLOTOP: 

PLOTPL: 

PLOTT: 

PLPLOT: 

PLRMS: 

PLO: 

PROFIL: 

PI: 

P2: 

Q: 

QOMU: 

QOMUP: 

QONU: 

QONUP: 


Number  of  ranges. 

Index  of  source  depth  of  profile  (may  be  input— see  Yarger  (1976)). 
Exponent  to  determine  spacing  of  wave  numbers  along  a path. 

2tr  frequency. 

1/3. 

Total  phase  change. 

Upper  phase  change. 

Lower  phase  change. 
tt  - 3.1415926535. 

Array  to  accumulate  real  vs.  range  part  of  property  loss  calculation 
for  coherent  phase. 

Array  to  accumulate  real  vs.  image  part  of  property  loss  calculation 
for  coherent  phase. 

Plot  option  (input — see  Yarger  (1976)). 

Plot  option  (input — see  Yarger  (1976)). 

Plot  option  (input — see  Yarger  (1976)). 

Plot  option  (input — see  Yarger  (1976)). 

Subroutine. 

Random  phase  propagation  loss  vs.  range. 

Minimum  dB  for  scale  of  property  loss  plot. 

Subroutine. 

Sign  of  receiver  depth  range  adjustment. 

Sign  of  source  depth  range  adjustment. 

Cycle  variable: 

Receiver  term  for  amplitude  calculation. 

Receiver  term  for  amplitude  calculation. 

Source  term  for  amplitude  calculation. 

Source  term  for  amplitude  calculation. 
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Minimum  range  (input — see  Yarger  (1976)). 


,1 

,1 

M 

! 

5 

.s 


RANGE: 


Range  array  in  yards. 


RAPLOT: 


Subroutine. 


RAYBLK: 


Common  block  for  plot  routines. 


Array  of  cycle  ranges  for  wave  numbers  along  a path. 


RCMAX: 


Maximum  RC  value  for  a single  J index,  all  4 paths. 


RCMIN:  Minimum  RC  value  for  a single  J index,  all  4 paths. 

RD:  Array  of  range  adjustments  to  cycle  range  for  receiver  depth. 


Linearly  Interpolated  RD  value. 

Real  part  of  a complex  term  for  each  mode. 


RESS: 


Surface  losses  in  dB. 


RMAX: 


Maximum  range  (input — see  Yarger  (1976)). 


Array  of  range  adjustments  to  cycle  range  for  source  depth. 


Linearly  Interpolated  RS  value. 


Subroutine. 


Surface  loss  values. 


Linear  Interpolated  R1  at  a particular  wave  number. 
Bottom  loss  array  dependent  on  wave  number. 


SALNTY: 


Salinity  in  0/00  for  XBT  conversion  (input — see  Yarger  (1976)). 


SETSMG: 


IGS  subroutine  to  set. 


Univac  trigonometric  function  for  sine;  argument  in  radians. 


SMLM: 


Smallest  mode  trapped. 


Univac  V function. 


SQUD: 


Subroutine. 


TESTB: 


Used  to  determine  range  cutoffs  for  "RAYMODE  Method". 


TESTQ: 


Used  to  determine  range  cutoffs  for  "SAYMODE  Method". 


TESTQP: 


Used  to  determine  range  cutoffs  for  "RAYMODE  Method". 
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THEDA: 

Array  of  angles  for  source  deviation  pattern  (input — see  Yarger 
(1976)). 

THEDA2: 

Array  of  angles  for  receiver  deviation  pattern  (input- 
(1976)). 

-see  Yarger 

THETA: 

Bottom  loss  angles  in  degrees  for  bottom  loss  table  (may  be  input — 
see  Yarger  (1976)). 

THETAS: 

Surface  loss  angles  in  degrees  for  surface  loss  table 
input — see  Yarger  (1976)). 

(may  be 

THREEH: 

Subroutine. 

TTPLOT: 

Subroutine. 

TWOPI: 

2tt  in  radians. 

WS: 

Wind  speed  (input — see  Yarger  (1976)). 

XIM: 

See  XRE  b>  low  except  for  imaginary  part. 

XRE: 

Term  containing  real  part  of  complex  expression  for  property  loss 
done  by  mode  summation. 

Z: 

Profile  depth  array  (input — see  Yarger  (1976)). 

ZB: 

Bottom  depth  (input — see  Yarger  (1976)). 

ZP: 

Depth  array  in  yards  of  adjusted  velocity  profile  used  internally. 

ZR: 

Receiver  depth  (input — see  Yarger  (1976)). 

ZRP: 

Receiver  depth  in  feet  or  meters  to  put  in  plot  key. 

ZS: 

Source  depth  (input — see  Yarger  (1976)). 

ZSP: 

Source  depth  in  feet  or  meters  to  put  in  plot  key. 

$ 
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• Submarine  Systems  Effectiveness  and 
Assessment  (SUBSEA) 

• Submarine  Active  Detection  Systems 
(SADS) 

• BQQ-5  Sonar  OMS 

• BQQ-6  Sonar  OMS 

• Sonar  In  Situ  Mode  Assessment  System 
(SIMAS) 

• Fleet  Mission  Library 

(U)  The  RAYMODE  X program  consists  of 
fifteen  subroutines  as  follows: 


£ 

I. 

RAMODX 

Propagation  Loss  Model. 
Main  Program 

& 

11. 

LEROY 

Temperature  to  Velocity 
Conversion 

» » 

i 

III. 

RSS 

Surface  Loss  Computation 

>5 

IV. 

MGSBL 

MGS  Bottom  Loss  Computa- 
tion 

V. 

THREEH 

Ray  Computations 

1 

VI. 

GETPRO 

Historical  Velocity  Pro- 
file Retrieval 

o 

VII. 

PROFIL 

Profile  Manipulation 

■ • 

VIII. 

LINTRP 

Linear  Interpolation 

5 

•V 

IX. 

BEAM 

Deviation  Loss  Computa- 
tion 

$: 

X. 

AMPTUD 

Amplitude  Calculation 

;s 

XI. 

SQUD 

Reduction  of  Trig  Func- 
tion Argument 

£ 

XII. 

DTAPLT 

Plot  SVP  and  Input  Loss 
Tables 

c 

XIII. 

RAPLOT 

Plot  Source  and/or  Re- 

ceiver versus  Range 

XIV. 

TTPLOT 

Plot  Travel  Time  versus 
Range 

N 

XV. 

PLPLOT 

Plot  Random  and/or  Coher- 

ent  Propagation  Loss  ver- 
sus Range 


(U)  In  the  user's  guide  of  Yarger  (1976) 
an  * marks  RAMODX  lines  of  code  and  plot- 
ting subroutines  DTAPLT,  RAPLOT,  TTPLOT , 
and  PLPLOT  to  be  removed  when  the  IGS 
graph  plotting  capability  is  not  avail- 
able. 

(U)  The  RAYMODE  model  has  error  stops 
generally  whenever  inputs  go  beyond  the 
allowable  range  as  given  in  Table  7-1  or 
exceed  dimension  parameters  established 
in  the  program.  If  an  error  exists  on 
the  HVP  (Historical  Velocity  Profile) 
tape  input,  an  error  message  will  print 
the  NTRAN  status.  Certain  error  condi- 
tions are  automatically  corrected  by  the 
program,  for  example,  a zero  profile 
gradient  or  Illegal  MGS  province. 

(U)  In  addition  to  accepting  a user- 
specified  sound  speed  profile,  RAYMODE  X 
at  NUSC/NL  can  access  historical  veloc- 
ity profile  (HVP)  data.  The  historical 
velocity  profile  data  on  magnetic  tape 
is  taken  from  NUSC  Technical  Documents 
5271,  5447,  5555,  and  6035  by  E.  Podes- 
zwa  that  contain  Sound  Speed  Profiles 
for  the  North  Pacific  Ocean,  North 
Atlantic  Ocean,  Indian  Ocean,  and  Nor- 
wegian Sea,  respectively. 

(U)  A chart  to  illustrate  the  path  and 
cycle  structure  as  well  as  signs  of  ray 
angles  in  RAYMODE  is  given  in  Figure 
2-1. 

3.0  (U)  The  Physics  of  the  RAYMODE  X Modal 
by  R.  Dnmport 

3.1  (U)  Introduction 

(U)  RAYMODE  X Is  a computer  program 
developed  by  G.  A.  Leiblger  at  the  New 
London  Laboratory  of  the  Naval  Under- 
water Systems  Center.  RAYMODE  Y calcu- 
lates the  acoustic  pressure  field  P(r,z) 
at  range  (r)  and  depth  (z)  due  to  a 
point  harmonic  source,  or  angular  fre- 
quency « , located  at  (r»o,  z-zg ) in  a 
p’ ecewise-layered  medium  (see  Fig.  3-1). 
TLansmlssion  loss  Is  then  calculated 
both  coherently  and  incoherently.  De- 
tails regarding  actual  running  of  the 
computer  program  can  be  found  in  Yarger 
(1976). 
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(U)  The  philosophy  adopted  by  Leiblger 
is  to  utilize  the  simpler  tools  of  ray 
theory  while  retaining  the  more  exact 
formulation  of  normal  mode  theory.  Two 
benefits  are  thus  realized:  (1)  it  is 

possible  to  interpret  mode  theory  ex- 
pressions in  a manner  similar  to  ray 
theory,  and  (2)  the  use  of  ray  theory 
simplifies  some  of  the  computational 
aspects  of  normal  mode  theory,  allowing 
for  considerable  savings  in  computer 
execution  time. 

(U)  The  basic  idea  of  RAYMODE  is  to  par- 
tition the  wavenumber  integral  solution 
for  the  acoustic  field  into  expressions 
which  individually  have  meaning  in  terms 
of  conventional  ray  theory.  This  is 
accomplished  by  using  the  velocity-depth 
profile  to  divide  the  wavenumber  domain 
into  regions  corresponding  to  surface 
duct  (SD),  convergence  zone  (CZ),  and 
bottom  bounce  (BB)  propagation  paths. 
Each  wavenumber  integral  is  then  expand- 
ed into  four  parts  corresponding  to  the 
four  rays  associated  with  the  upgoing/ 
downgoing  eigenrays  at  the  source/re- 
ceiver. These  upgoing/downgoing  ray  in- 
tegrals are  then  numerically  evaluated 
by  either  normal  mode  theory  or  via  a 
multipath  expansion.  Therefore,  each 
propagation  path  consists  of  four  parts 
and  the  total  field  is  the  result  of 
summing  all  paths  (i.e.,  SD  + CZ  + BB). 

(U)  This  synopsis  attempts  to  present 
the  theory  and  approximations  leading  up 
to  the  equations  which  are  ultitn  ttely 
used  in  the  RAYMODE  X program. 

3.2  (U)  Theory 

Integral  Solution  (U) 

(U)  The  acoustic  pressure  P(r,z)  is 
found  by  applying  the  Fourier-Bessel 
transform  to  the  Helmholtz  wave  equation 
in  cylindrical  coordinates  (r,  9 , z). 


where  azin  ial  symmetry  is  assumed  and 
c(z)  is  the  sound  speed  as  a function  of 
depth.  In  addition,  a time  factor  of  exp 
(iu»t)  has  been  used.  The  Fourier- 
Bessel  solution  of  (3-1)  is  given  by 


Mr,*,*  ) 


Iw 


/nu,v 


t)  h(  (cotdt. 


0-2) 


o 

where  Hq  is  the  Hankel  function  of  order 
zero  and  G is  the  depth  dependent 
Green's  function.  In  general  G is  a very 
complicated  expression  for  a piecewise 
continuous  medium.  However,  if  it  is  as- 
sumed that  sound  speed  profile  disconti- 
nuities do  not  backscatter  appreciable 
energy,  then  the  Green's  function  can  be 
written  as 


[*(*,)  + K*  f<*,)Kf(*>  + **  g(*>]  °"3) 


and  f and  g are  linearly  independent 
traveling  wave  solutions  of  the  separat- 
ed wave  equation  in  the  depth  variable 
z. 


(3-4) 


where 


q(*> 


z 

and  the  wave  phase  4>  N is  defined  by 
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(U)  The  depth  zjj  is  the  phase  refer- 
ence depth  for  the  wave  functions  in 
zg,  and  is  either  the  ocean  surface  or 
an  upper  turning  point  depth  correspond- 
ing to  a downward  refracted  wave*  Like- 
wise, Zty  is  the  phase  reference  depth 
for  the  wave  functions  in  z and  is 
either  the  bottom  of  a surface  duct,  the 
ocean  bottom  depth  or  a lower  turning 
point  depth  corresponding  to  an  upward 
refracted  wave.  The  superscripts  K and  N 
on  the  reflection  coefficients  Ruand  R^j, 
refer  to  the  above  reference  depths, 
Zfc  and  Zty,  respectively.  Thus,  R^  is 
evaluated  at  z^  and  R^j  is  evaluated  at 
zjj.  The  subscripts  u and  d on  the  re- 
flection coefficients  r£  and  R^  refer  to 
the  direction  (up  or  down)  of  propaga- 
tion of  the  wave  functions  that  are  re- 
flected at  Zfc  and  z^.  Thus  R^f(z) 
represents  the  reflected  wave  of  the 
upward  traveling  wave  g(z).  Likewise 
R^g(z)  is  the  reflected  wave  of  the 
downward  traveling  wave  f(z).  The  re- 
flection coefficients  must  therefore  be 
chosen  such  that  the  boundary  conditions 
at  Zfc  and  zjj  are  satisfied. 

, Depth  Dependent  Solutions  for  a Segmented 
Velocity  Profile  (U) 

(U)  The  RAYMODE  method  assumes  that  the 
sound  speed  profile  c(z)  can  be  approxi- 
mated by  segments  such  that  q(z)  is  a 
linear  function  of  z in  each  segment 
(layer).  Therefore,  within  each  layer 
segment,  independent  traveling  wave  so- 
lutions of  (3-4)  are  exactly  given  by 
Hankel  functions  of  order  1/3  and  may  be 
written  in  WKB  form  as 


depth  for  which  q(z)»0)  within  any 
layer,  so  that  V(z)  is  bounded  near  to 
and  at  the  turning  point. 

Continuity  Conditions  and  Reflection 
Coefficients  (U) 

(U)  Due  to  the  layering  of  the  profile, 
the  solutions  f(z)  and  g(z),  and  their 
derivatives,  should  be  continuous  across 
all  interfaces.  This  matching  of  the  so- 
lutions (for  exponential  layers)  is  done 
exactly  in  the  Fast  Field  Program  (FFP). 
However,  in  RAYMODE  X this  matching  is 
only  done  for  the  interface  at  the  bot- 
tom of  a surface  duct,  where  R^  is  re- 
placed by  a reflection  coefficient  ob- 
tained from  the  continuity  conditions 
for  a bilinear  profile.  Everywhere  else 
Leibiger  makes  the  approximation  that  a 
wave  is  totally  transmitted  across  all 
interfaces  until  either  a layer  is 
reached  where  a turning  point  exists 
the  wave  interacts  with  either  the  ocean 
surface  or  bottom.  At  a turning  point 
the  reflection  coefficient  is  assigned 
unit  magnitude  and  a tt/2  phase  shift. 
At  the  ocean  surface  a plane  wave  re- 
flection coefficient  is  assumed  with  a 
— ti  phase  shift  and  a magnitude  obtained 
from  a semiempirical  formula  based  on 
the  works  of  Beckman  and  Spizzichino 
(i5S3)  and  Marsh  and  Schulkin  (1962). 
Surface  loss  (i.e.,  -20  log^oiR|)  ob- 

tained from  this  magnitude  consists  of 
two  parts:  a high  frequency  loss,  SLj, 
and  a low  frequency  loss,  SL2*  The 
surface  loss  SLj  is  given  by 


f(z)  • V ( z) axp 


SLj  ■ -20  Lo*l0  (1-V3r 


where 


z 

“V  (i) exp  H q'hlidl  | 


where  the  amplitude  V(z)  and  its  complex 
conjugate  V*(z)  are  given  in  Appendix 
3A.  It  is  important  to  note  that  the 
profile  is  partitioned  so  that  there  is, 
at  most,  only  one  turning  point  (i.e., 
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larger  than  .99,  then  V3  is  set  equal  to 
.99.  The  loss  SL2  is  given  by 


SL  - -20  log 


10 


.3  + • 


(1  + .01  (2 


WS  . 10 


V 


| (3-10) 


(U)  The  above  assumptions  regarding  the 
reflection  coefficients  are  equivalent 
to  approximating  solutions  of  equation 
(3-4)  over  [z^,  Zjj]  by  generalized 
WKB  solutions  which  are  valid  at  turning 
points. 


where  f is  the  frequency  in  Hz.  Details 
regarding  the  derivation  of  (3-8),  (3-9) 
and  (3-10)  are  given  in  Appendix  IIIB. 
(A  user  of  RAYMODE  X can  also  specify 
his  own  surface  loss.) 

(U)  At  the  ocean  bottom  the  phase  of  the 
reflection  coefficient  is  assumed  zero 
and  the  magnitude  is  obtained  from  a set 
of  modified  MGS  curves  or  by  directly 
entering  one’s  own  bottom  loss  table.  If 
nothing  is  specified,  zero  dB  loss  is 
assumed. 


Integration  Limits  (U) 

(U)  For  computational  purposes  the  inte- 
gration limits  on  equation  (3-2)  are  ap- 
proximated and  the  integration  performed 
piecewise.  For  example,  for  the  profile 
given  by  Figure  3-2,  with  source  and 
receiver  both  in  the  surface  duct,  the 
integration  limit n are  partitioned  90  as 
to  correspond  to  ray  angles  appropriate 
for  surface  duct  (SD)  paths,  convergence 
zone  (CZ)  paths  and  bottom  bounce  (BB) 
paths.  Thus, 


(U)  Due  to  the  fact  that  the  MGS  curves 
were  originally  meant  for  frequencies 
greater  than  1000  Hz  (specifically  3500 
Hz),  the  RAYMODE  X program  modifies 
these  curves  for  frequencies  less  than 
1000  Hz.  In  particular  for  100  Hz  and 
below,  a single  curve  derived  by  Christ- 
ensen, et  al.  (1973)  is  used  for  all  MGS 
provinces.  This  curve  yields  bottom  loss 
BLj  (9)  as  a function  of  grazing  angle 
(9)  in  degress  and  is  given  by 

BL| (8)  - -3.11  + .4040  - 4.98  x lO'V  + 2.89  x 10'V  (3-11) 

- 9.0  x 10"V 


jXo  0 -*2 

/-(  -{  */ 

— tj  £l  e3 

(3-13) 

(SD)  (CZ)  (BB) 

or 

pap  + p 4-  p 

SD  CZ  BB 

where 

- 

W/Ci  (1-0, 1,2) 

and 

the 

velocitites 

Cq,  Cj.  ad  C2 

are 

de- 

fined 

by  the 

profile  in  Fig. 

3-2. 

5 3 

is  defined  by  the  largest  source  (graz- 
ing) angle  9 s to  be  considered  by  us- 
ing the  relation, 


(U)  For  frequencies  between  100  Hz  and 
1000  Hz  bottom  Iose  BL2  (9)  as  a func- 
tion of  frequency  (f),  province  (P)  and 
grazing  angle  is  found  by  Interpolating 
between  the  100  Hz  curve  and  the  MGS 
curve  for  the  particular  province. 

(U)  Explicitly,  BL2  (0)  is  given  by 

PL2(6)  - BLt  (0)  + [bL3(6  ,P)  - BLj  (0)]  X Log^f.Olf)  <3-'2) 


where  BL3  (6,P)  represents  the  appro- 
priate MGS  curve  as  a function  of 
grazing  angle  9 and  province  £. 


5j  -(—* V086.  (3-14) 

\ Source  / 

(U)  In  general  the  limits  of  integration 
are  approximated  as  follows:  start  with 

the  Min  jc0,  csource  creceiver  | = 
cMin  an^  search  for  the  next  velocity 
maximum  cj_;  this  defines  the  first  set 
of  limits,  ( ?i,  Cmsx)*  where  5Max« 
w/c^in*  From  cj  to  the  next  velocity 
maximum  yields  the  secon ! integra- 
tion interval  (52,  51).  This  process  is 
continued  until  the  last  velocity  maxi- 
mum is  reached,  which  is  usually  the 
water  velocity  eg  at  the  ocean  bottom. 
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The  final  set  of  integration  limits  are 
then  (£  final*  ^b)  where  final  is 
determined  by  equation  (3-14)  when  the 
largest  desired  source  angle  is  speci- 
fied . 

(U)  From  the  above  description  of  the 
integration  limits  It  Is  important  to 
note  that,  for  surface  duct  situations, 
complex  source  angles  are  included  in 
the  integration,  thus  allowing  for 
leakage  (diffraction)  effects. 

jilpgoing  and  Downgoing  Paths  (U) 

(U)  Each  of  the  £ -partitioned  integrals 
for  the  pressure  P is  now  rewritten  as 
the  sura  of  four  integrals 

(P  = P +p  + p +p  ) 

' BB  BB!  BB2  BB3  BB4' 

obtained  when  the  numerator  of  the 
Green's  function  (3-3)  is  expanded  in 
(3-2).  These  four  integrals  can  be 
identified  with  the  four  ray  paths 
associated  with  the  upgoing/downgoiog 
eigenrays  at  the  source  and  receiver, 
for  that  particular  5 -partition.  For 
example,  when  the  £ -partition  is  for 
bottom  bounce  paths,  the  four  integrals 
correspond  to  the  paths  shown  in  Fig. 
3-3.  The  easiest  way  to  see  this  associ- 
ation is  to  use  first  order  WKB  forms 
for  the  solutions  * and  g,  and  then 
integrate  each  integral  by  stationary 
phase.  The  stationary  phase  points  are 
simply  the  target  (eigen-)  rays  for  the 
four  upgoing/downgoing  rays  at  the 
8ource/rec.eiver  combination.  This  is 
interesting,  but  only  gives  a way  to 
obtain  ray  theory  from  a wave  theory. 
Leibiger  has  generalized  the  above 
picture  by  using  generalized  WKB  forms 
for  f(z)  and  g(z)  and  the  integrating 
the  above  ray  path  integrals  very 
accurately. 

Numerical  Evaluation  of  Integrals  (U) 

(U)  Each  of  the  four  integrals  is 
evaluated  by  either  normal  mode  theory 
(when  the  number  of  modes  within  the 
^-partition  is  small)  or  via  a multipath 
expansion.  Leibiger  has  found  that  when 


the  number  of  modes  exceed  ten  it  is 
computationally  expedient  to  calculate 
the  field  integrals  by  the  multipath  ex- 
pansion method.  This  allows  for  a con- 
siderable savings  in  computer  execution 
time  since  the  multipath  expansion 
method  effectively  sums  the  higher  order 
modes  in  one  fast  integration.  In  both 
instances  the  integrals  (for  same 
^-partition,  ^b)»  to  be  evalu- 
ated are  of  the  form, 


,,  -/■  2S* 


tie'1  <*i-*k‘  + 5r> 


(1  -RkR« 
u d 


A(z-V°  ' '1\1Gr  ) V <S>  V<*> 


where  for  illustrative  purposes  only 
of  the  above  four  integrals  is 
considered.  Also  the  first  terra  in  the 
asymptotic  expansion  of  H^)((r)  is 
assumed. 

Normal  Mode  Evaluation  (U) 

(U)  When  normal  mode  theory  is  utilized, 
the  singularities  K ^ associated  with 
the  modes  are  assumed  to  be  simple  poles 
obtained  by  solving  the  equation 


V -('  - R"Itde'21^N)-  0 


(U)  The  normal  mode  residue  associated 
with  (3-15)  then  becomes 


A(z, -if**  - *ZfJ  + £ r) 

P(  ) = 2»i  L Jwf-  e k k m ' 'X1) 

1 - 


Details  regarding  the  numerical  deter- 
mination of  the  complex  eigenvalues 
and  specific  evaluation  of  the  normali- 
zation term  -j^j^cau  be  found  in  Appen- 
dix 3C. 
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Multipath  Expansion  (U)  ) 

(U)  For  the  multipath  approach,  the 
denominator  of  (3-15)  is  expanded  so 
that 

Eg 

P(  >1  ° j?ot  A<***,!0(R' “ ^)J  '*k‘"2;1*kN  + 5r]  O'18) 


The  interval  of  integration  is  divided 
into  a number  of  unequal  sections  based 
upon  the  number  of  rays  traced  (which  is 
an  input  parameter).  The  value  of  £ for 
some  of  these  sub-sections  is  suffi- 
ciently far  from  a stationary  phase 
point  so  that  their  contribution  is  ex- 
cluded. Although  stationary  phase  tech- 
niques are  not  used,  the  stationary 
phase  points  are  available  from  ray  cal- 
culations. The  justification  for  neg- 
lecting such  sub-sections  is  involved 
with  the  spiral-like  nature  of  the  cumu- 
lative result  for  the  field.  When  a sub- 
section does  not  meet  this  criteria  the 
phase  term  is  approximated  by  a quadrat- 
ic expression  in  5.  The  amplitude  of  the 
kernel  is  assumed  to  be  slowly  varying 
over  the  sub-interval  so  that  it  can  be 
evaluated  at  an  interior  point  and  re- 
moved outside  the  integral.  The  result- 
ing integral  can  be  expressed  in  terms 
of  Fresnel  integrals  by  a suitable 
transformation.  The  Fresnel  integrals 
are  then  evaluated  numerically.  Fresnel 
integrals  can  be  used  to  evaluate  the 
integrals  over  the  sub-intervals  because 
of  the  smallness  of  the  integration  in- 
terval. If  larger  intervals  were  used 
the  computation  time  would  increase  due 
to  the  need  of  evaluating  incomplete 
Airy  functions. 

(U)  The  series  given  by  (3-18)  repre- 
sents the  multipath  expansion  of  (3-15). 
The  advantage  of  (3-18)  is  that  not  only 
are  the  upgoing-downgoing  paths  deline- 
ated but  the  number  of  cycles  that  a ray 
(wave)  undergoes  is  counted  by  the  index 
j.  In  RAYMODE  X the  number  of  cycles 
necessary  for  a ray  path  to  reach  a 
given  range  is  calculated  so  that  the 
infinite  series  in  (3-18)  may  be  approx- 
imated by  sudsing  only  a few  j's.  For 


example,  if  the  5-partition  corresponds 
to  convergence  zone  paths  and  one  is  ,i- 
terested  only  in  the  first  CZ,  ' -..in 
terms  for  j»l  would  only  yield  *i'e 
structure  effects  ra  a propagation  loss 
versus  range  curve. 

Transmission  Loss  (U) 

(U)  The  acoustic  pressure  field  as  de- 
termined from  (3-17)  and  (3-18)  is  modi- 
fied by  a beampattern  attenuation  factor 
characterizing  the  off-axis  beam  posi- 
tion of  an  equivalent  ray.  This  beampat- 
tern attenuation  is  similarly  applied  to 
the  other  three  pressure  components 

P(  )2>  P(  )3  and  P(  >4  for 
each  5 -partition.  The  real  and  imaginary 
components  of  pressure, 

Re  |P(  )i|  and  Im  jP  ( )±\ 

are  then  used  to  form  a rms  intensity 
and  a coherent  intensity  for  each  range 
point  desired  and  for  each  ^-partition . 

(U)  Transmission  loss  (relative  to  unit 
intensity  at  unit  distance  from  source) 
is  then  calculated  from  the  incoherent 
and  coherent  intensity  sums.  For  exam- 
ple, in  the  case  described  in  Fig.  3-2, 
the  transmission  loss  TL  is  given  by 

l ' * V 

TL  (Coherent)  - -10  lo*10  + M^l  + )2 

+(i4tiB,PsDi| + in,,w + 


TL  - (incoherent)  - -10  l°g10  j ^ 53  [CRe ji>SD  |)2  + (Im|PSD  j)2 

+ (Re.]  P^])2  * (ImlP^l)2 

+ <R*ipBBl|>a+ (im,x,)2]! 

+ or  (3-20) 

where  is  Thorp's  attenuation 

coefficient . 
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3.3  (U)  Summery  of  Buie  Assumption 

• Velocity  profile  fit  with  segments 
such  that  the  index  of  refraction 
squared  is  a linear  function  of  depth. 

5 • Multiple  reflections  due  to  velocity 

* discontinuities  ignored  except  in  sur- 
face duct  situations. 

| • Plane  wave  reflection  coefficients  as- 

sumed. 

| • Surface  a d bottom  loss  expressions 

1 assume  the  validity  of  experimental 

^ data,  which  may  be  questionable. 

) • Only  a finite  number  of  ray  cycles  are 

considered  in  the  multipath  evaluation 
| of  the  pressure  integrals. 

• Velocity  profile  does  not  change  with 

| range . 

* • Harmonic  source  assumed. 


(TJ)  In  this  case  RAYMODE  does  not  prop- 
erly account  for  the  surfact  duct  (SD) 
contribution  because  only  trapped  modes 
are  presently  used  in  the  UNIVAC  1108 
version.  Trapped  modes  are  here  defined 
to  mean  those  modes  whose  eigenvalues 
are  real  and  lie  between  £ l * w/ci 
and  £2  " w/c2*  where  c^  is  the 

velocity  at  the  layer  depth  and  C2  is 
the  larger  of  the  source/receiver 
velocities.  For  the  lower  frequencies 
one  must  also  allow  for  the  leaky  modes, 
i.e.,  modes  whose  eigenvalues  are  com- 
plex. The  leaky  modes  correspond  to 
eigenvalues  between  £ 0 - w/c0  and 

5*  - (“> /c source)  coa  9s>  where  co 

is  the  velocity  at  the  surface  and  08 
is  the  largest  source  (grazing)  angle. 
The  SD  part  is  described  as  allowing  Lot 
both  trapped  and  leaky  modes.  Thai  's, 
the  SD  integration  interval  extends  from 
£0  to  £»  . However,  this  version  does 
not  exist  on  the  UNIVAC  1108  although  it 
does  exist  on  the  HP9845  and  the 
Tektronix  4051  computers. 
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• Density  of  1 assumed. 

• Only  one  source  and  receiver  allowed 
per  run. 

• Constant  bottom  depth. 

3.4  (U)  Suggested  Test  Coses  for  RAYMODE  X 

(U)  The  RAYMODE  X program  has  been  com- 
pared against  both  experimental  data  and 
other  computer  programs.  In  general  the 
comparisons  tend  to  validate  most  of  the 
assumptions  listed  in  Sections  3.3  of 
this  synopsis.  However,  there  are  a few 
cases  that  may  possibly  cause  RAYMODE  X 
to  give  unsatisfactory  answers.  Results 
for  two  cases  are  given  below: 

(U)  Case  I:  Cross  layer  surface  duct 

problem  with  source  in  the  duct  at  250 
ft  and  receiver  below  the  duct  at  450  ft 
for  frequencies  of  10,  50  and  100  Hz. 
The  velocity  profile  for  this  case  is 
shown  in  Figure  3-4.  Bottom  loss  is 
given  by  FN0C  type  5. 


(U)  For  the  in-layer  case  (source  At 
76.2  m,  receiver  at  45.72  m)  coherent 
RAYMODE  (Figs.  3-5  to  3-7)  predicts  no 
trapped  SD  modes  at  50  Hz  and  100  Hz. 
Therefore,  only  the  direct  path  and  bot- 
tom bounce  contribute  at  these  frequen- 
cies. In  order  to  evaluate  RAYMODE,  Fast 
Field  Program  (FFP)  predictions  were 
made  (Figs.  3-8  to  3-10)  since  the  FFP 
considers  all  modes.  When  RAYMODE  is 
compared  with  the  FFP,  one  observes  that 
in  the  fifst  bottom  bounce  region  the 
leaky  mode  contribution  is  masked  by  the 
bottom  bounce  energy.  However,  beyond 
the  first  bottom  bounce  region  the  leaky 
modes  are  the  significant  contribution. 

(U)  For  the  in-layer  case  at  200  Hz 
RAYMODE  predicts  that  one  (unattenuated) 
trapped  mode  exists.  When  compared  with 
the  FFP  it  is  clear  that  RAYMODE  pre- 
dicts much  more  trapping  (approximately 
10  dB)  than  actually  exists  because  even 
at  200  Hz  the  true  trapped  mode  has  a 
significant  imaginary  component.  This 
imaginary  part  attenuates  the  SD  mode 
contribution.  (Note:  Incoherent  RAYMODE 
results  are  given  in  Figures  3-11  to 
3-13.) 
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Figure  3-4.  (U)  Velocity 
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(U)  In  the  ctv-js-layer  case  (source  at 
76-i  m,  receiver  at  13/.  16  m)  RAYMC.  E 
(Figs.  3-14  to  3-16)  showed  basically 
good  results  compared  to  the  FFP  (Figs. 
3-17  to  3-19)  for  all  frequencies  con- 
sidered despite  the  fact  that  it  does 
not  allow  for  below  layer  leakage 
effects.  That  is,  in  the  present  1108 
version  it  does  not  matter  in  the  cross 
layer  situation  whether  or  not  there  are 
trapped  modes  in  the  duct  since  no  ener- 
gy is  allowed  to  escape.  (Note:  RAYilODE 
incoherent  results  are  given  in  Figures 
3-20  to  3-22.) 

(U)  In  summary,  Case  I shows  that  the 
present  UNIVAC  1108  version  of  RAYMODE 
can  yield  poor  surface  duct  results  at 
the  lower  frequencies.  This  deficiency 
in  RAYMODE  has  been  recognized  by  Lei- 
biger  for  some  time,  and  has  been  cor- 
rected in  the  HP  and  Tektronix  versions. 

(U)  Case  II:  Depressed  sound  channel 
problem  with  source  and  receiver  both  in 
the  depressed  channel  at  depths  of  270 
ft  and  220  ft,  respectively,  and  for 
frequencies  of  10,  30,  100  and  300  Hz. 
Velocity  profile  for  this  case  is  shown 
in  Figure  3-23.  Bottom  loss  is  given  in 
Table  3-1. 


Table  3-1.  (U)  Bottom  Loss  Versus  Graz- 
ing Angle  for  Test  Case  II 


Angle  (degrees) 
0 

Bottom  Loss  (dB) 
BL(O) 

0.0 

0.0 

0.5 

2.8 

1.0 

5.6 

1.5 

8.6 

2.0 

11.5 

2.5 

14.2 

3.0 

16.7 

3.5 

18.8 

4.0 

20.5 

5.0 

23.0 

6.0 

24.4 

8.0 

25.7 

10.0 

26.1 

15.0 

26.4 

90.0 

26.4 
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(U)  Case  II  is  concerned  with  low  fre- 
quency propagation  in  an  environment 
where  the  velocity  profile  (Fig.  3-23) 
exhibits  a depressed  channel  above  the 
S0FAR  channel.  The  source  and  receiver 
are  located  near  the  axis  of  the 
depressed  channel.  For  a frequency  of  10 
Kz,  RAYMODE  (Fig.  3-24)  doer  not  see  the 
depressed  channel  (depressed  modes  below 
phase-integral  cutoff)  but  only  seme 
combination  of  convergence  zone  (CZ)  and 
bottom  bounce  energy.  The  FFP  (Fig. 
3-27),  however,  seems  to  indicate 
trapping  within  the  depressed  channel. 
When  the  frequency  is  increased  to  100 
Hz,  RAYMODE  (Fig.*  3-25)  still  does  not 
calculate  depressed  modes  but  predicts 
more  loss  than  at  10  Hz,  which  is  diffi- 
cult to  understand  since  there  should  be 
more  trapping  at  100  Hz.  The  FFP  for  100 
Hz  (Fig.  3-28)  does  predict  more  trap- 
ping within  the  shallow  channel  as  well 
as  some  CZ  energy  at  40  km  and  80  km.  At 
300  Hz  RAYMODE  (Fig.  3-26)  appears  now 
to  be  predicting  one  totally  trapped 
depressed  mode  along  with  some  CZ  energy 
at  around  43  km  and  86  km.  The  FFP  for 
300  Hz  (Fig.  3-29)  predicts  strong  de- 
pressed channel  trapping  (approximately 
12  dB  less  loss  than  RAYMODE)  with  only 
hints  of  the  CZ  energy  at  40  km  and  80 
km.  (Note:  RAYMODE  incoherent  results 
are  given  in  Figs.  3~30  to  3-32.)  In 
summery  RAYMODE  does  not  properly 
account  for  depressed  channal  propaga- 
tion at  the  lower  frequencies  (<300  Hz) . 
The  exact  reasons  are  not  clear  but  seem 
to  be  related  to  the  fact  that  in  RAY- 
MODE there  is  no  consideration  of 
partial  trapping  of  energy  within  the 
depressed  channel > This  low  frequency 
trapping,  however,  is  masked  many  times 
by  low  loss  bottom  bounce  energy. 

Appendix  3A.  Depth  Dependent  Solutions 
for  a Segmented  Velocity  Profile  (U) 

(U)  The  RAYMODE  method  assumes  that  the 
sound  speed  profile  c(z)  can  be  approxi- 
mated by  segments  such  that  q(z)  is  a 
linear  function  of  z in  each  segment 
(layer).  Therefore,  within  each  layer 
segment,  independent  traveling  wave 
solutions  of  (3.4  are)  exactly  given  by 
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Figure  3-24.  (U)  Case  II.  RAYMODE  Coherent.  Source  Depth  = 82  m.  Receiver  Depth 

67  m.  Frequency  = 10  Hz. 
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Figure  3-32.  (U)  Case  II.  RAYMODE  Incoherent.  Source  Depth  = 82  m.  Receiver  Depth 

67  m.  Frequency  = 300  Hz. 
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where  HW3  are  Hankel  functions  of 
order  1/3;  T and  V(z)  are  complex 
quantities  defined  by 


Appendix  3B.  Surface  Scattering  Loss  (U) 

(U)  Surface  loss  in  RAYMODE  X assumes 
two  scattering  mechanisms:  a high  fre- 
quency, large  roughness  loss,  SL^,  and 
a low  frequency  loss,  SL2.  The  total 
surface  scattering  loss  is  then  given  by 
the  sum  SL}  + SL2.  It  is  not  clear 
what  the  physical  basis  is  for  these  two 
processes.  1 suggest  that  it  would  be 
profitable  if  this  aspect  of  RAYMODE  X 
be  examined  (especially  since  a user  of 
RAYMODE  X may  input  his  own  surface  loss 
table).  The  large  roughness  SLi  is 
given  by 


n 

SL,  - -10  lo*ln  1-  I / oco.*rd*rdOr  (3B- 


■ t (4 ) 


T*  and  V*(Z)  represent  the  complex  con- 
jugates of  T and  V(Z) . The  phase  is 
defined  as 


which  is  identical  to  (3-5)  but  for 
variable  reference  depth. 

(U)  The  traveling  wave  solutions  f and  g 
can  thus  be  written  in  generalized  WKB 
form,  , 


f(x)  - V(i)  exp  I- 


j-1 qV*)  dr| 


g(i)  - V 


z 

*(z)  exp|  ij 


and  still  be  exact  within  any  given 
layer. 


where  the  angles  are  defined  by  the  ge- 
ometry shown  in  Figure  3B-1.  Equation 
(3B-1)  is  based  on  the  assumption  that 
an  intensity  scattering  (reflection) 
coefficient  R should  obey  the  equation 


}"rn 

■ I I ocoa^d^dfi^  - 


(U)  Therefore,  if  surface  loss  is  de- 
fined by  -10  logjolRl»  then  (3B-2) 
yields  the  loss  formula  given  by  (3B-1). 
The  scattering  coefficient  o is  the 
usual  one  given  when  the  Fraunhofer 
phase  approximation  is  assumed  (for  ex- 
ample, see  Beckman  and  Splzzichino, 
(1963).  In  the  large  roughness  limit  o can 
be  written  in  terms  of  the  distribution 
of  slopes  as 


-cot2e0f(*r,er>  cot28n fz 

* (ain>,  + iind^)2 


where  the  mean  square  slope  is  h tan2/30 
and 

cot20o  - [2( .003  + 2.6  x 10~3WS)]"1 
A coa2$.  - 2 coa$,cos$  co#8  + cos2$ 

t(W  ■ (.In*,  + «in*r)2  ( 


(ain$^  + Bin$f) 
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and  WS  is  the  windspeed  in  knots.  When 
(3B-3)  is  used  in  (3B-1)  and  the  inte- 
grations performed,  the  surface  loss 
Slq  is  approximately  given  by 

SL  » -20  log  (1-V3)11  (3B-5) 


However,  there  is  something  unusual 
about  (3B-8),  because  according  to 
Sverdrup  and  Munk  (1947)  the  signifi- 
cant wave  height  H3/3  is  given  by 

H1/3(ft)  - 2.32  x 10“2WS2 (knot) , (3B-10) 


where 


and  if  the  Pierson-Moskowitz  (1964) 
spectrum  is  used,  H1/3  becomes 


V3  “ maximum  of 


oln9- 

alnd 

2 


faln6 

8 


(3B-6) 


and  a ■ cot  30»  6 is  the  grazing  angle  in 
radians.  If  V3  is  larger  than  .99,  then 
V3  is  set  equal  to  <99. 


H1/3(ft)  - 1.86  x 10'2WS2(knot).  (3B-11) 

(U)  The  usual  relationship  between  the 
average  wave  height  H and  the  signifi- 
cant wave  height  H^/3  is  given  by 

H - ,625H1/3.  (3B-12) 


(U)  The  "low  frequency"  surface  loss 
term  SL2  is  based  on  surface  duct 
scattering  data  given  by  Marsh  and 
Schulkin  (1962).  A suitable  equation 
describing  this  experimental  data  is 
given  by 


(U)  When  (3B-12)  is  used  in  (3B-10)  and 
(3B-11)  one  obtains 

H(ft)  - 1.45  X 10“2WS2 (knot)  (3B-13) 

(Sverdrup-Munk) 


SL 


2 


-20  log 


10 


and 

(3B-7) 


H(ft)  - 1.16  x 10_2WS2 (knot)  (3B-14) 

(Plerson-Motkovltz) 


where  f represents  the  frequency  in  Hz 
and  H is  the  average  wave  height  in 
feet.  Eugene  Podeszwa*  at  NUSC/NL  has 
analyzed  some  results  of  Vine  and  Volk- 
mann  (Woods  Hole,  1950)  and  arrived  at 
the  following  relationship  between  wave 
height  (H)  in  feet  and  wind  speed,  WS, 
in  knots: 

H - 2.04  x 10~2WS2  . (3B-8) 

(U)  When  this  formula  of  Podeszwa’s  is 
used  in  (3B-7),  the  surface  loss  SL2 
in  terms  of  wind  speed  becomes 
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(U)  If  the  average  height  given  by  OB- 
13)  and  (3B-14)  is  compared  with 
Podeszwa' s (3B-8),  it  would  appear  that 
Podeszwa 's  "wave  height"  corresponds 
more  closely  to  the  significant  wave 
height  H1/3. 

(U)  However,  in  describing  the  Marsh- 
Schulkin  data  it  is  necessary  to  use  the 
average  wave  height.  Because  of  the 
above  discrepancy,  it  is  felt  that  the 
use  of  an  alternate  formula  should  be 
considered . 

Appendix  3C.  Normal  Mode  Eigenvalues  and 
Eigenfunction  Normalization  (U) 

(U)  The  complex  eigenvalues  £m  associ- 
ated with  the  normal  mode  (residue)  ex- 
pansion (3-17)  are  found  from  the  zeros 
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of  the  denominator  of  (3-15).  Explicitly 
they  are  assumed  to  be  single  poles  ob- 
tained by  solving  the  equation 


' (‘  - 


(3C-1) 


If  the  reflection  coefficients  are 
defined  in  terms  of  magnitude  and  phase 
(0K  and  9^)  then  the  above  eigen- 
value equation  becomes 


W M - .2U,<N) 


9,  + e„  + 2 Re(t'N)  - 2mtr,  a - 0,1,2. 
It  N k 


(3C-2) 


The  real  part  of  the  eigenvalue , t ra  - 
Re  j5m}  is  found  by  solving  (3C-2) 
without  recourse  to  iterative  methods 
which  results  in  a considerable  savings 
In  execution  time.  To  accomplish  this 
the  extreme  mode  numbers  (%,Ng)  are 
determined  by  substituting  5 a an<*  € B 
(the  end  points  of  a given  5 -partition) 
into  the  second  part  of  (3C-2)  and  solv- 
ing form-  In  addition,  differentiation 
yields 


dm 

di 


1_ 
2 n 


Rc  (m) 


!!n! 

« J 


(3C-3) 


where 


Rc(m)  - -i-^2  S.u]N))|, 


(3C-4) 


is  the  cycle  range^  associated  with  the 
mode  eigenvalue  5 m.  The  value  of 
Rc(m)  is  available  in  closed  form  due 
to  the  assumed  sound  speed  variation 
within  each  layer.  If  one  assumes  that 
the  phase  of  the  reflection  coefficients 
is  a slowly  varying  function  of  5,  the 
above  equation  may  then  also  be  solved 
for  the  extreme  values  5 a and  C g* 
Thus  a curve  of  5 m vs.  mode  number  ra 
passing  through  the  end  points  ( 5 a» 
Na)>  ( 5 B»  nB>  and  having  the 

slopes  prescribed  above,  can  be  obtained. 


This  curve  is  approximated  by  a cubic 
polynomial  and  the  unknown  coefficients 
are  determined  by  interpolation.  This 
expression  is  then  evaluated  for  each 
integer  such  that  NA<5m<Ni$  for  a direct 

(non-iterative)  evaluation  of  the  real 
part  of  the  eigenvalues. 

(U)  The  imaginary  part  of  the  eigen- 
values 5m  is  approximately  given  by 


ImU  ) »■ 


. Rk 
-lnAs_u_ 


Rc(m) 


(3C-5) 


The  normal  mode  residue  expansion  asso- 
ciated with  (3-15)  then  becomes 


B A(z,Z  }C  ) j / j2  A*8  >.  r r\ 

P.  . (r,z,z  ) = 2wi  £ -g-  ‘ *k  + 5mr) 

( ) , “ 


-r,  3w  | 

m-NA  ac 


(3C-6) 


The  normalization  term  in  the  denomina- 
tor can  be  approximately  written  as 


Rr  (m) 


OC-7) 


and  further  approximated  by  retaining 
only  the  term  involving  R^.  The  final 
form  of  the  residue  normal  mode  series 
is  obtained  by  expanding  the  phase  terms 
(e.g.,  <j)£  and  <j£s)  ) in  a Taylor  series 

about  5m  and  retaining  only  the  first 
two  terms,  since  for  trapped  modes 
(Im|5m|  is  small.  Then 


r(  ) 


(r,z,z  ) * 2 


k Tk  m 


_ ^ (a (t, z ; £ ) , 

. I * a -i ( 

/ ^ I Rc(m)  e 
m-N  u 
A 

x eIm(F-mt  Cr(zk>z)  - r(VV  ' r]J 


(3C-8) 
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4.0  (U)  Running  Tims 


(U)  The  running  times  for  the  cases 
examined  in  the  RAYMODE  X evaluation  are 
given  in  Table  4-1.  Running  times  are 
the  total  of  simultaneously  exercising 
the  incoherent  and  coherent  phase  addi- 
tion options.  All  running  times  were 


obtained  on  the  UNIVAC  1108  computer  at 
the  New  London  Laboratory.  Input  varia- 
bles for  these  runs  were  allowed  to  as- 
sume default  variables  (except  for  most 
environmental  inputs. 


Table  4-1.  (U)  RAYMODE  X Run  Times  on  Univac  1108 


Datn  Set 

Case 

Source 
Depth  (m) 

Receiver 
Depth  (m) 

Frequency 

(Hz) 

Number  of 
Points 

Run  Time 
(Sec) 

HAYS-MURPHY 

I 

24.4 

137.2 

35.0 

400 

9.6 

II 

24.4 

137.2 

67.5 

400 

51.9 

III 

24.4 

137.2 

100.0 

400 

53.7 

IV 

24.4 

200.0 

400 

56.0 

V 

24.4 

SS 

35.0 

400 

10.3 

VI 

24.4 

106. 1 

100.0 

400 

55.2 

F/.SOR 

I (FIG) 

6.1 

37.0 

150U.0 

200 

13.4 

23.0 

37.0 

1500.0 

200 

3.1 

III  (THORN; 

23.0 

37.0 

1500.0 

200 

4.8 

IV  (REDWOOD) 

6.1 

37.0 

1500.0 

200 

14.5 

PARKA 

I 

152.4 

91.4 

50.0 

400 

29.2 

II 

152.4 

91 .4 

400.0 

300 

54.1 

BEARING  STAKE 

I 

91.0 

496.0 

25.0 

300 

6.5 

II 

91.0 

1685.0 

25.0 

300 

6.9 

III 

91.0 

3320.0 

25.0 

300 

6.6 

IV 

91.0 

3350.0 

25.0 

300 

S9 

. 

V 

13.0 

496.0 

1400 

300 

S3 

VI 

18.0 

1685.0 

140.0 

300 

5.1 

VII 

18.0 

3320.0 

140.0 

300 

5.4 

18.0 

3350.0 

140.0 

300 

5.3 

18.0 

496.0 

290.0 

300 

5.6 

WSKSKM 

18.0 

1685.0 

290.0 

300 

XI 

18.0 

3320.0 

290.0 

XII 

18.0 

3350.0 

290.0 

5.4 

JOAoT 

I 

mm 

18-3 

3700.0 

200 

16.5 

II 

wmm 

79.2 

3700.0 

18.1 



III 

mam 

163.1 

3700.0 

200 

16.7 

IV 

mi 

18.3 

3700.0 

200 

23.6 

V 

■sS 

79.2 

3700.0 

200 

24.6 
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Data  Set 


JOAST  (cont.) 


Gulf  of  Alas 


I (Run  140) 


II  (Run  140) 


III  (Run  143) 


IV  (Run  143) 


V (Run  124) 


VI  (Run  124) 


VII  (Run  112A) 


III  (Run  112A) 


IX  (Run  112B) 


30.5 


30.5 


30.5 


30.5 


30.5 


30.5 


1066.8 


1066.8 


1066.8 


Receiver 

Frequency 

Number  of 

Run  Time 

Depth  (m) 

(Hz) 

Points 

(Sec) 

163.1 

3700.0 

200 

24.2 

18.3 

3700.0 

200 

15.8 

79.2 

3700.0 

200 

15.4 

163.1 

3700.0 

200 

15.9 

18.3 

3700.0 

200 

15.8 

163.1 

3700.0 

200 

15.8 

18.3 

3700.0 

200 

16.2 

79.2  J 

3700.0 

200 

16.2 

304.8 

3700.0 

200 

16.2 

17.0 

400.0 

400 

21.5 

112.0 

400.0 

400 

22.3 

43.0 

1000.0 

400 

22.3 

112.0 

1000.0 

400 

22.4 

6.0 

1500.0 

400 

21.9 

59.0 

1500.0 

400 

25.9 

6.0 

2500.0 

400 

41.1 

59.0 

2500.0 

400 

17.0 

3500.0 

400 

1 

112.0 

3500.0 

400 

22.4 

17.0 

5000.0 

400 

17.5 

112.0 

5000.0 

400 

18.2 

30.5 

1500.0 

250 

30.3 

304.8 

1500.0 

250 

32.1 

30.5 

1500.0 

250 

31.3 

304.8 

1500.0 

250 

31.5 

30.5 

1500.0 

250 

30.0 

304.8 

1500.0 

250 

32.4 

30.5 

2500.0 

250 

28.4 

304.8 

2500.0 

250 

26.2 

30.5 

2500.0 

250 

27.2 
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Data  Set 

Case 

— 

Source 
Depth  '-t) 

Receiver 
Depth  (m) 

Frequency 

(Hz) 

Number  of 
Points 

Run  Time 
(Sec) 

Gulf  of  Alaska 
(cont.) 

X (Run  1123) 

1066.8 

304.8 

2500.0 

250 

28.1 

XI  (Run  107) 

304.8 

30.5 

2500.0 

2ZC 

18.0 

XII  (Run  107) 

304.8 

304.8 

2500.0 

250 

17.7 

XIII  (Run  108) 

304.8 

30.5 

2500.0 

250 

26.8 

XIV  (Run  108) 

304.8 

304.8 

2500.0 

250 

28.9 

LORAD 

IA 

15.2 

30.5 

530.0 

300 

23.0 

IB 

15.2 

30.5 

530.0 

300 

23.0 

I.C 

15.2 

30.5 

530.0 

300 

24.9 

ID 

15.2 

30.5 

530.0 

300 

24.9 

IE 

15.2 

30.5 

530.0 

300 

24.3 

IF 

15.2 

30.5 

530.0 

300 

24.3 

IG 

15.2 

30.5 

530.0 

300 

24.3 

IIA 

15.2 

304.8 

530.0 

300 

22.4 

I IB 

15.2 

304.8 

530.0 

300 

22.4 

IIC 

15.2 

304.8 

530.0 

300 

24.7 

IIC 

15.2 

304.8 

530.  r 

- 

300 

24.7 

HE 

15.2 

304.8 

530.. 

300 

24.2 

IIF 

15.2 

304.8 

530.0 

300 

24.2 

IIG 

15.2 

304.8 

530.0 

300 

24.2 
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5.0  (U)  Core  Storage  Requirements 

(U)  The  number  of  decimal  words  of  core 
storage  for  RAYMODE  is  presented  in 
Table  5-1  in  four  ways  for  each  of  two 
different  sets  of  parameters  used  in 
dimensioning  data  arrays.  Core  storage 
was  cclculated  both  with  and  without  the 
1GS  plotting  routines,  but  including 
other  normal  input/output  functions. 
These  numbers  have  been  broken  down  into 


data  and  code  size  for  the  RAYMODE  pro- 
gram routines  alone  and  for  the  program 
with  the  UNIVAC  system  routines  neces- 
sary to  execute  it.  The  program  stored 
is  sized  according  to  the  first  set  of 
parameters  to  accommodate  even  large 
problems;  the  second  set  shows  the  ef- 
fect of  reducing  these  parameters  to  a 
still  reasonable  size. 


Table  5-1.  (U)  RAYMODE  Stortge  Size 


Method 

Parameters 

Code 

Data 

Total 

With  plots  and  systems 

dimensioned  for 

22213 

11927 

34140 

With  plots,  without  system 

400  ranges 

7084 

5547 

12631 

Without  plots,  with  system 

and  50  modes  and 

10236 

7083 

17319 

Without  plots  and  system 

ray  pts 

3928 

4533 

8461 

With  plots  and  system 

dimensioned  for 

22213 

10502 

32715 

With  plots,  without  system 

200  ranges 

7084 

4122 

11206 

Without  plot 8 , with  system 

and  25  modes  and 

10236 

5658 

15894 

Without  plot 8 and  system 

ray  pts 

3928 

3108 

7036 
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6,0  (U)  Program  How 


(U)  The  following  flowchart  for  the 
RAYMODE  model  was  obtained  from  D.  F. 
Yarger  of  the  Naval  Underwater  Systems 
Center,  New  London  Laboratory. 


(1  Of  5) 
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7.0  (U)  RAYMODE  X Input! 

(U)  RAYMODE  X inputs  and  their  discus- 
sion by  Yarger  (1976)  are  presented  be- 
low. Some  editorial  liberties,  including 
rearrangement  and  rewriting  of  some  text 
as  necessitated  by  the  format  of  this 
report,  have  been  undertaken.  Also,  de- 
tailed discussion  of  inputting  histori- 
cal velocity  profiles  (HVPs)  has  been 
omitted . 

7.1  (U)  RAY  HI  ODE  Control  Card  and  Data  Deck 
Requirements 


the  EXEC  VIII  executive  system.  The 
propagation  model  is  stored  on  cata- 
logued file  TEN*RAYM0D  and  should  be 
executed  from  the  absolute  element  RAMO 
DXA.  The  control  stream  for  running  RAY- 
MODE is  shown  in  Figure  7-1,  which  also 
contains  a sample  multi-case  data  deck. 

(U)  If  a historical  velocity  profile 
(HVP)  is  to  be  called,  the  user  must  as- 
sign the  correct  tape-to-tape  unit  1 be- 
fore execution.  The  tape  number  to  be 


(U)  RAYMODE  is  presently  running  on  the 
UNIVAC  1108  computer  at  NUSC/NL  under 


@RUN 

@ASG,A 

@ASG,T 

@USE 

@XQT 

$ INPUTS 


SEND 

SINPUTS 

SEND 


. . . , run-time, max-pages 
TEN*RAYMOD 

HVTAPE,U,U79«  (or  U544) 

1..HVTAPE. 

TEN*RAYMOD . RAMODXA 
NORTH  ATLANTIC  - WINTER 


(run-time  .5n,  max-pages  20n  for  n cases) 


lonly  if  HVP 
>to  be  used 


UN  .'VAC 
EXEC  VIII 
Control 
Cards 


(heading  card) 
(in  col.  2-8) 


I OCEAN  - 2,  IPROFL  - 10,  ISEASN  - 4, 

ZS  - 20. , ZR  - 200. , 

F - 500., 

MGSOP  - 3,  WS  - 15., 

R-  1000.,  DELTAR  - 1000.,  RMAX  - 100000., 
NORTH  ATLANTIC  - WINTER 


(historical  profile) 
(source  & receiver) 
(frequency) 

(bottom  loss  A wind  speed) 
(range  determinants) 

(in  col.  2-5) 


) 


Sample 
Case  1 


F - 1000., 


(change  frequency) 


Case 


2 


NORTH  ATLANTIC  - WINTER 


$ INPUTS 

• . 

SEND 

0FIN 


UNCLASSIFIED 


Case 


n 


A 


B 


Figure  7-1.  (U)  RAYMODE  Control  and  Data  Deck  Structure 
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referenced  depends  on  the  unit  system  of 
the  data,  English  or  metric;  where 

Tape  U769  contains  HVP's  in  ft  and 
ft/s  (English  system). 

Tape  U544  contains  HVP’s  in  m and 
m/s  (metric  system). 

The  tape  will  automatically  position  it- 
self to  the  proper  HVP  even  when  calling 
several  different  ones  in  the  same 
execution. 

(U)  Following  the  UNIVAC  control  in- 
structions as  outlined  in  Figure  7-1, 
each  separate  case  requires 

• a heading  card  in  format  12A6 

• the  required  NAMELIST  inputs  inserted 
between  the  indicators  $INPUTS  and 
SEND. 

(U)  The  input  information  to  be  supplied 
by  the  user  consists  of: 

(1)  Sound  velocity  profile:  input 

directly  or  form  historical  velocity 
profile  tape. 

(2)  Source  and  Receiver  depths:  ZS, 
ZR  or  indices  NU,  MU. 

(3)  Frequency:  F. 

(4)  Horizontal  range  determinants: 
R,  RJAX,  DELTAR 

(5)  Bottom  loss:  KGS  province  number 
of  fable. 

(6)  Wind  speed:  WS. 

(7)  Source  and/or  Receiver  beam  pat- 
terns (defaulted  omnidirectional). 

C8)  Program  controls  (defaulted  for 
general  use). 

(9)  Output  options  (defaulted 
"ou").* 

♦Editors  Note:  See  section  7.2  for 

RAYMODE  X output  options. 


(U)  A discussion  of  the  use  of  NAMELIST 
and  the  particular  RAYMODE  inputs  within 
each  of  the  above  categories  is  included 
in  this  section.  The  NAMELIST  method 
offers  an  advantage  when  running  more 
than  one  case.  For  the  first  case,  the 
user  needs  to  include  only  those  input 
variables  not  supplied  by  defaults  plus 
the  ones  for  which  the  defaults  are  not 
suitable  for  his  application;  for  each 
case  after  the  first,  the  user  needs  to 
include  only  those  input  variables  which 
differ  from  the  ones  previously  set.  The 
illustration  in  Figure  7-1  shows  that  a 
typical  case  in  the  English  unit  system 
(default)  using  any  of  the  historical 
velocity  profiles  (3  inputs)  for  any 
sc.  rce/receiver  combination  (2),  fre- 
quency (1),  MGS  province  (1),  wind  speed 
(1),  and  set  of  ranges  (3)  with  omnidi- 
rectional source  and  receiver  (default) 
and  calling  for  all  available  output 
(default)  requires  exactly  11  inputs  to 
be  defined  with  the  appropriate  numeri- 
cal values.  To  change  only  frequency  for 
a second  case  merely  requires  another 
heading  card  and  another  NAMELIST 
$INPUTS  deck  with  the  new  frequency  as- 
signed. Therefore,  although  the  NAMELIST 
list  of  inputs  may  appear  somewhat  for- 
midable because  it  includes  up  to  three 
varieble  names  (eight  of  which  are 
arrays  with  up  to  50  entries  each),  the 
usual  cases  are  handled  quite  simply 
with  the  program  defaults.  The  extensive 
list  of  Inputs,  however,  allows  for 
great  flexibility  of  the  program  control 
in  special  applications.  User  inputs  are 
not  modified  internally.  A summery  of 
RAYMODE  inputs  with  their  limits  and  de- 
faults is  provided  in  Table  7-1. 

(U)  The  specific  input  requirements  for 
execution  of  each  RAYMODE  case  are  dis- 
cussed below: 

1.(U)  Heading  Card 

(U)  The  first  card  required  for  each 
RAYMODE  data  case  is  an  alphanumeric 
heading  card  read  in  format  12A6.  The 
comment  on  this  card  is  intended  to  be 
the  title  of  the  case  and  will  be 
printed  at  the  beginning  of  the  printed 
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Table  7-1  (cont.)*  (U)  Summary  of  Raymode  Inputs 
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output  and  at  the  bottom  of  any  selected 
plots.  For  best  results  on  the  plots, 
any  comment  should  be  centered  within 
the  first  72  columns.  If  no  title  Is 
desired,  a blank  card  must  be  used. 

2.  (U)  NAMELIST  $iNPUTS  Deck 

(U)  A detailed  discussion  of  the  NAME- 
LIST  procedure  Is  contained  In  the  FOR- 
TRAN V manual,  reference  (e).  In  brief, 
the  first  card  of  the  NAMELIST  deck  will 
be  $TNPUTS  In  columns  2-8  and  the  last 
will  be  $FND  in  columns  2-5;  in  be- 
tween, the  Inputs  are  defined  In  equa- 
tion form  and  separated  by  commas,  for 
example : 

(a)  in  the  case  of  a constant 

F-50., 

and 

(b)  In  the  case  of  an  array  using  a 
subscript  and  listing  the 
elements 

THETA(l)  - 0.,  5.,  10.,  ...»  55.,  60., 

(U)  The  inputs  defined  in  this  way  af- 
ford convenience  to  the  user  in  checking 
data  for  later  review.  For  this  reason  a 
print  of  the  data  deck  itself  is  mean- 
ingful, hence  is  provided  as  part  of  the 
output  as  shown  in  the  examples  below. 
Variables  within  the  NAMELIST  deck  may 
be  assigned  in  any  order  and  need  not  be 
In  special  fields  or  on  separate  cards 
with  the  exception  that  the  use  of  col- 
umn 1 should  be  avoided.  Each  veriable 
must  be  of  the  appropriate  type,  integer 
or  real..  For  RAYMODE  each  input  variable 
id  typed  according  to  the  name  rule 
(i.e.,  real  unless  beginning  with  the 
letters  I-N)  with  the  exception  of  the  4 
plot  opt  Loos  which  are  integers.  Arrays 
are  distinguished  by  the  subscript  (1). 

(U)  The  particular  input  categories  are 
discussed  in  the  following  section. 


(U)  A sound  velocity  profile  may  be 
input  directly  from  cards  listing  such 
point  from  surface  to  bottom  using  N,  Z, 
C,  SALNTY  oir  be  accessed  from  a histori- 
cal velocity  profile  tape  by  providing 
the  three  parameter  values  I0CEAN, 
IPR0FL,  ISEASN. 

a.  (U)  Selection  of  doits:  METRIC 

(U)  The  systems  of  units  to  be  used  for 
input  and  output,  is  controlled  by  the 
integer  METRIC. 

METRIC  ■ 0 implies  the  English  system 
(default) 

* 1 implies  the  Metric  system. 

b.  (U)  Directly  Input  Profile:  N,  Z,  C, 

SALNTY 

(U)  When  using  the  English  system  if  the 
depth  units  are  feet,  a profile  may  be 
made  up  of  a mixture  of  velocities  in 
ft/s  and  temperatures  in  °F  in  order  to 
easily  merge  bathythermograph  (6T)  data 
with  an  SVP.  When  using  the  metric 
system,  a profile  may  be  made  up  of  a 
mixture  of  velocities  in  m/s  and  temp- 
eratures in  °C  for  the  same  reason.  The 
program  automatically  converts  tempera- 
ture data  to  sound  velocity  measure  in 
the  appropriate  units  to  complete  the 
SVP.  Adjacent  equal  velocities  or  temp- 
eratures are  accepted  by  the  program  but 
are  slightly  modified  for  internal  use 
to  avoid  a zero  gradient  condition.  The 
user  inputs  for  a directly  inserted  SVP 
are: 

N - no.  cf  pts  in  SVP,  2<N<47. 

Z(l)  - depth  entries  from  surface  to 
bottom  in  increasing  order  in 
ft  or  yd  if  METRIC  = 0,  in  m 
if  METRIC  - 1. 

C(l)  * velocity  entries  corresponding 
to  Z depths  in  ft/s  or  yd/ s if 
METRIC  - 0,  in  m/s  if  METkJC  - 
1,  c>£  temperatures  in  °F  if 

depth  in  ft  and  METRIC  * 0, 
in  °C  if  METRIC  - 1,  C±  < 100. 
when  temperature. 
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SALNTY  ■ salinity  in  parts  per  thousand 
for  conversion  of  temperatures 
to  velocities  if  any  tempera- 
tures are  entered  in  C array 
above  (default  35). 

3.  (U)  Access  Historical  SVP:  IOCEAN, 

IPROFL,  ISEASN 

(U)  The  historical  velocity  profile  tape 
was  developed  from  the  seasonal  profiles 
selected  by  E.  Podeszwa  (February,  1976; 
April,  1976). 

IOCEAN  ■ ocean  code  (default  0). 

_<  0 implies  that  SVP  is  input 
directly  as  in  b above. 

> 0 Implies  that  historical  pro- 
file is  to  be  used. 

* 1 selects  the  North  Pacific 
Ocean 

* 1 selects  the  North  Atlantic 
Ocean. 

IPROFL  ■ profile  index  a3  obtained  from 
area  charts. 

ISEASN  - season  index  1 <ISEASN  <4. 

- 1 for  winter  (Jan-Mar) . 

■ 2 for  spring  (Apr-Jun). 

■ 3 for  summer  (Jul-Sep). 

» 4 for  fall  (Oct-Dec). 

4.  (U)  Bottom  Depth  ZB 

(U)  Bottom  depth  may  be  inserted  direct- 
ly in  the  same  units  as  the  profile 
depths  input.  This  input  is  particularly 
required  for  use  with  a historical 
velocity  profile,  since  a fixed  bottom 
depth  of  21,000  feet  (or  6400  meters)  is 
assumed  for  the  stored  profiles.  The 
program  truncates  or  extends  the  stored 
HVP  to  conform  to  the  selected  bottom 
depth.  Since  charted  depths  are  often 
uncorrected,  the  user  is  responsible  for 
any  correction  to  bottom  depth  before 
eucering  it.  The  input  bottom  depth  will 
be  used  as  the  last  point  of  the  profile 
when  >0.;  if  >0.,  the  last  profile  depth 
Is  assumed  to  be  the  bottom. 


ZB  - bottom  depth  in  ft,  yd,  or  m cor- 
responding to  depth  units. 

£ 0.  implies  Z(N)  is  bottom  depth. 

> 0.  implies  bottom  depth  input  spe- 
cifically. 

Source  and  Receiver  Depths  ( U ) 

(U)  Source  and  receiver  depths  may  be 
identified  by  the  indices  of  these 
depths  in  the  profile  using  NO  or  MU  or 
by  the  actual  depths  using  ZS  and  ZR, 
respectively. 

a.  (U)  Indices  of  Source  and  Receiver 
Depths:  NU,  MU 

(U)  When  the  source  and  receiver  depths 
are  both  points  in  the  profile,  the  user 
may  identify  these  depths  by  using  the 
appropriate  indices. 

NU  ■ index  of  source  depth,  1>NU>N. 

■ 0 implies  ZS  and  ZR  to  be  used. 

MU  ■ index  of  receiver  depth,  1>MU>N. 

b.  (U)  Source  and  Receiver  Depths:  ZS,ZR 

(U)  When  either  source  or  receiver 
depths  are  not  points  in  the  profile  or 
when  using  a historical  profile,  the 
user  may  insert  the  depths  directly  in 
the  same  units  as  the  profile  depths  in- 
put. When  using  this  method,  set  NU  * 0. 

ZS  ■ source  depth  in  ft,  yd,  or  m corre- 
sponding to  depth  units, 
0.<ZS<bottom  depth. 

ZR  » receiver  depth  in  same  units  as  ZS, 
0.<ZR<bottom  depth - 

Frequency:  F (U) 

(U)  RAYMODE  has  been  run  for  frequencies 
as  low  as  10  Hz  and  as  high  as  100  kHz. 

F ■ frequency  in  Hz,  F>0. 
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Range  Determinants:  R,  RMAX,  DELTAR  (U) 

(U)  The  program  will  accommodate  up  to 
400  ranges. 

R * range  minimum  in  yd  if  h.  '’RIC-O. 

in  m if  METRIC-1, R>0. 

RMAX  - range  maximum  in  same  units  as 
R,RMAX>R. 

DELTAR  - range  increment  in  sam*  units 
as  R,  DELTAR>0. 

Bottom  Loss  ( U ) 

(U)  MGS  Bottom  Loss  Province  MGS0P 

(U)  When  an  MGS  bottom  loss  province 
number  is  entered,  the  program  generates 
a bottom  loss  table  for  grazing  angles 
evary  2n  from  0°  to  the  maximum  sonar 
angle  ANGLE  discussed  later. 

MGS0P  - MGS  bottom  loss  province 

O<MGS0P<9  (default  0). 

b.  (U)  Bottom  Loss  Table;  ITAB , THETA, BL 

(U)  Set  MGS0P  « 0 when  an  input  bottom 
loss  table  is  desired. 

, TAB  - no.  of  pts  in  inpu^  bottom  loss 
table,  0<TTAB<50  (default  0). 

THETA  (1)  - angles  in  degrees  for  bottom 
loss  table  to  be  supplied 
from  0"  t"  ANGLE  in  increas- 
ing order. 

BL(1)  - bottom  losses  in  dB>0.  corre- 
sponding to  THETA  angles. 

Wind  Speed:  WS  (U) 

(u)  For  a positive  wind  speed  input,  the 
program  will  generate  a surface  loss  vs. 
grazing  angle  table  every  2°  from  0°  to 
ANGLE. 

WS  - wind  speed  in  knots,  WSX>. 

(default  0). 


Beam  Patterns  for  Source  and/or  Receiver  (U) 

(U)  Vertical  beam  patterns  la  the  form 
of  deviation  loss  vs.  angle  may  be 
applied  to  the  source  and/or  receiver. 
The  absence  of  a deviation  loss  table 
indicates  an  omnidirectional  pattern. 
In  defining  the  deviation  loss  table 
downward-directed  angles  from  the  hori- 
zontal (0°)  are  interpreted  by  the  pro- 
gram as  being  positive,  upward-directed 
angles  as  negative. 

a.  (U)  Source  Deviation  Loss  Table:  IDL, 
THEDA,  DL 

IDL  - no.  of  pts  in  source  deviation 
table,  0<IDL<50  (default  0). 

THEDA(l)  » angles  in  degrees  for  source 
deviation  loss  table  from 
-ANGLE  to  ANGLE  in  increasing 
order . 

DL(1)  * source  deviation  losses  in  dB>0. 
corresponding  to  THEDA  angles. 

b.  (U)  Receiver  Deviation  Loss  Table: 
JDL,  THEDA2,  DL2 

JDL  »•  no.  of  pts  in  receiver  deviation 
table,  0<JDT,<50  (default  0). 

DL2(1)  - receiver  deviation  losses  in 
dB>0.  corresponding  ro  THEDA2 
angles . 

Program  Control r '!) 

(U)  As  previouc^y  noted,  che  default 
values  of  the  program  control  parameters 
and  those  controls  which  are  internally 
generated  by  the  program  itself  are 
adequate  for  usual  applications  of  the 
model.  Program  control  is  provided  for 
the  user  who  is  sufficiently  familiar 
with  certain  pertinent  aspects  of  the 
model  to  employ  successfully  the  input 
parameters  described  in  this  section. 
These  controls  general!  v vary  the  amount 
of  computation  performed  in  the  execu- 
tion of  the  model  and  therefore  impact 
the  computation  time  and  associated 
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coat.  Some  of  the  controls  affect  the  imum,  LAMMIN<LAMDA<50  (default 

accuracy  of  the  computation  so  that  the  -1). 

user  may  run  the  model  in  accordance 


with  his  own  specific  accuracy/running 
time  requirements.  Generally,  the  great- 
er the  required  accuracy  the  greater  the 
running  time.  Thus,  if  a user  has  many 
cases  to  run  and  is  not  overly  concerned 
with  providing  the  maximum  attainable 
accuracy,  he  may  set  the  program  con- 
trols to  process  each  case  rapidly, 
thereby  keeping  the  total  accumulated 
run-time  for  all  cases  within  some  ac- 
ceptable limit. 

(U)  Another  use  of  the  program  controls 
involves  the  possibility  of  separate 
examination  of  distinct  contributors  to 
the  total  acoustic  field.  For  example, 
the  user  may  desire  to  compare  the  rela- 
tive effect  of  surface  4ict  and  conver- 
gence zone  propagation  on  the  total 
propagation  loss,  or  to  compare  the 
relative  contributions  of  single  normal 
modes  or  subsets  of  modes.  The  employ- 
ment of  program  controls  in  the  manner 
described  below  allows  such  special 
model  applications  to  be  processed. 

a.  (U)  Control  of  Number  of  Ray  Cycles: 
LAMMIN,  LAMDA,  LAMDAB 

(U)  The  RAYMODE  program  will  compute  the 
minimum  and  maximum  number  of  ray  cycles 
necessary  for  both  nonbottom  bounce 
(e.g.,  surface  duct  or  convergence  zone) 
and  bottom  bounce  propagation  loss  com- 
putations by  using  selected  cycle  range 
values.  The  user  may  wish  to  assign  his 
own  cycle  values  to  reduce  execution 
time  or  to  isolate  selected  propagation 
paths,  e.g.,  first  convergence  zone). 
Negative  cycle  inputs  indicate  the 
cycles  are  under  program  control  only. 

LAMMIN  » minimum  number  of  ray  cycles 
for  all  propagation  modes  if 
computed  minimum,  LAMMIN<50 
(default  -1). 

LAMDA  » maximum  number  of  ray  cycles  for 
all  propagation  modes  other  than 
bottom  bounce  if  > computed  raax- 


LAMDAB  ■ maximum  number  of  ray  cycles 
for  bottom  bounce  if  < computed 
maximum,  LAMMIN<LAMDAB<50  (de- 
fault -1). 

(U)  Therefore,  if  the  user  enters  LAMMIN 
■ 1,  LAMDA  " 1,  LAMDAB  ■ 2,  the  program 
will  compute  one  ray  cycle  for  all  non- 
bottombounce  propagation  paths  with  one 
and  two  bottom  reflections  only.  Setting 
LAMMIN  ■ 0 will  cause  the  program  to 
include  the  short  range  direct  and 
surface-reflected  propagation  paths 
(cycle  0). 

b.  (U)  Control  of  Source  Angle  Limits: 
ANGLO,  ANGLE' 

(U)  Selection  of  angle  limits  will  re- 
strict the  computation  of  propagation 
loss  to  that  portion  of  the  source  ener- 
gy emitted  within  the  ray  angle  limits 
defined  by  ANGLO,  ANGLE.  Thus,  the  user 
can  isolate  that  portion  of  the  total 
propagation  structure  influencing  sur- 
face duct  (SD),  convergence  zone  (CZ), 
and  bottom  bounce  (BB)  for  separate  ex- 
amination. 

(U)  For  user  convenience  ANGLO  and/or 
ANGLE  may  be  entered  in  terms  of  sound 
velocities  in  the  same  units  as  the  pro- 
file velocity  inputs.  The  chart  below 
illustrates  how  to  select  these  limits 
for  particular  combinations  of 
propagation  paths. 

Cq  ■ max  sound  velocity  on  the  SVP  be- 
tween and  including  ZS  and  ZR. 

Cl  ■ surface  layer  velocity. 

Cjj  = bottom  velocity. 

Cj[  * maximum  velocity  on  entire  SVP. 

ANGLO  =*  positive  and  negative  minimum 
sonar  angle  in  degrees, 
0.<ANGL0<90.  (default  0). 
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or  velocity  limit  in  same  units 
~ as  profile  velocities. 

ANGLE  - positive  and  negative  maximum 
aonar  angle  in  degrees, 
ANGL0<ANGLE<90 . (default  60). 

c.  (U)  Control  of  Modes  Processed: 
MINM0D,  MAXM0D,  MO 

(U)  The  specification  of  inputs  to  con- 
trol mode  selection  will  result  in  the 
program  evaluation  of  propagation  loss 
for  only  the  normal  modes  between  and 
including  the  indicated  limit  values. 
This  control  can  therefore  be  used  to 
examine  the  field  strength  of  individual 
modes  or  subsets  of  the  total  mode  set 
by  setting  limits  MINM0D  and  MAXM0D 
when  the  total  number  of  modes  trapped 
is  <MO,  the  cutoff  for  normal  mode  proc- 
essing. A value  of  0 for  MAXM0D  will 
result  in  the  program  computing  relative 
mode  numbers  from  MINM0D  to  the  larg- 
est mode  number  for  trapped  modes 
(NMODE ) . 


(U)  Thus,  if  MINM0D  - 5,  MAXM0D  = 5 and 
MO  - 10,  the  program  will  treat  only  the 
fifth  mode  of  the  sequence  1,  . .., 

NMODE,  ignoring  the  remaining  NMODE  - 1 
modes  in  the  determination  of  propaga- 
tion loss. 

MINM0D  * first  relative  normal  mode 
processed  by  mode  summation, 
KMINM05O  (default  1). 

MAXM0D  ■ last  relative  normal  mode 
processed  by  mode  summation, 
MAXM0D  <50  (default  0). 

M0  ■ maximum  number  of  modes  p-o^essed 
by  mode  summation,  1<M0<50  (default 
10). 

d.  (U)  Number  of  Points  in  Ray  Tables: 
NL 

(U)  Selection  ->f  the  number  of  points  in 
the  ray  tables  controlled  by  the  integer 
input  NL  greatly  influences  running  time 
and,  to  a lesser  extent,  the  accuracy  of 
the  model  results.  A preliminary  step 
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leading  to  propagation  loss  evaluation 
is  the  construction  of  tabular  data 
involving  repeated  application  of  a ray 
tracing  algorithm.  Once  obtained,  the 
information  is  processed  sequentially  as 
required  in  the  propagation  loss  rou- 
tine. The  larger  the  NL  the  more 
accurate  the  propagation  loss  results 
but,  as  expected,  the  longer  the  run- 
time. Values  of  NL  as  high  as  50  pro- 
vide the  most  precise  results  but  are 
rarely  needed.  Accuracy  suffers  little 
in  reducing  NL  to  20  while  improving 
run-time  by  more  than  a factor  of  two. 
For  typical  applications  NL  ■ 10  is 
almost  always  satisfactory,  with  errors 
in  propagation  loss  being  limited  to 
approximately  0.5  dB;  hence,  this  value 
of  NL  is  used  as  the  default  condition. 
Reducing  NL  to  values  of  five  or  less 
yields  exceedingly  rapid  execution  times 
with  errors  expected  in  the  1-2  dB 
region  at  worst. 

NL  - number  of  points  in  ray  tables, 
2<NL<50  (default  10). 

(U)  An  example  of  the  input  data  deck 
for  a RAYMODE  run  using  a Pacific  Ocean 
sound  8 peed  profile  from  'cards  in  the 
English  system  is  given  in  Table  7-2. 
Table  7-3  gives  the  input  data  deck  for 
an  example  which  accesses  a historical 
sound  8 peed  profile  from  the  North 
Atlantic  Ocean  in  the  metric  system. 

7.2  (U)  RAYMODE  X Outputs 

(U)  RAYMODE  X outputs  and  their  discus- 
sion as  presented  by  Yarger  (1976)  are 
presented  below.  Some  editorial  liber- 
ties, including  rearrangement  and  re- 
writing of  some  text  as  necessitated  by 
the  format  of  this  report,  have  been  un- 
dertaken. 

(U)  Tables  7-4  and  7-5  give  the  printed 
outputs  for  the  examples  corresponding 
to  the  inputs  of  Tables  7-2  and  7-3,  re- 
spectively. The  RAYMODE  model  will  pro- 
vide printout  of  the  selected  inputs  and 
tabulated  propagation  loss  versus  range 
and  optionally  provide  ray  information 
tables.  The  printout  will  be  titled  with 


the  comments  used  on  the  heading  card. 
The  velocity  profile  input  from  cards  or 
accessed  from  tape  will  be  printed  in 
its  original  units  followed  by  source, 
receiver,  and  bottom  depths;  then  the 
velocity  profile  used  internally  by  the 
program  in  yd  with  source,  receiver, 
and  bottom  depths  inserted  is  printed. 
The  other  prints  related  to  input  are 
frequency,  wind  speed  with  the  generated 
surface  loss  table  if  WS>0.;  MGS  prov- 
ince and  its  generated  bottom  loss  table 
or  the  bottom  loss  table  input  directly, 
if  any;  the  source  deviation  loss  table, 
if  any;  the  receiver  deviation  loss 
table,  if  any;  sonar  angle  limits;  mode 
inputs;  and,  finally,  a computed  refer- 
ence velocity  Cq  used  for  ray  computa- 
tions. Printouts  of  the  ray  information 
will  follow  if  this  option  is  chosen. 
Finally  the  output  table  of  Range  (R)  in 
kyd  or  kilometers,  depending  on  the  unit 
system  used,  versus  coherent  phase  prop- 
agation loss  (PL),  which  represents  the 
coherent  summation  of  all  computed  prop- 
agation paths  and  random  phase  propaga- 
tion loss  (PLRMS) , which  represents  an 
intensity  summation  over  all  propagation 
paths  in  units  of  decibels  relative  to 
one  yard  from  the  source  (dB//l  yd),  is 
produced  for  the  number  of  ranges  se- 
lected. A maximum  of  twenty  printed 
pages  will  be  produced  per  case. 

(U)  The  model  will  also  automatically 
plot  propagation  loss  versus  range  when 
used  in  conjunction  with  the  NUSC/NL 
graph  plotting  facility,  the  Information 
International  FR  80,  which  utilizes  the 
Integrated  Graphics  System  (IGS)  soft- 
ware; there  are  four  other  available 
plots  that  the  user  may  select.  The  exe- 
cution time  addition  for  the  total  plot 
set  is  less  than  ten  seconds.  The  set  of 
available  output  plots  furnishes  the 
most  convenient  means  for  Intepretation 
of  RAYMODE  results;  however,  the  poten- 
tial user  employing  the  model  at  a com- 
puting facility  without  the  IGS  graph 
plotting  capability  must  rely  on  the 
printed  version  of  the  output.  Such  a 
user  will  have  little  difficulty  in  in- 
terpreting the  printed  output  after  some 
scrutinization  of  the  sample  plots  in 
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The  input  data  deck  tor  Example  1 is  presented  below: 


PARKA  Data 


I 


& 


i 


P 


$INPUTS 

N = 29, 

Z(  1)  = 0., 50. ,100. ,131. 2, 164. ,180. 5, 200. ,229. 7, 246. 1,262. 5, 
278.9,300. ,400. ,500. ,600. ,800. ,1200. ,1312.3,1968.5,2500. , 
3280.8,4921.2,6561.7,8202.1,9842.5,11482.9,13123.3,16404.2, 
18044.6, 

C(l)  = 5022.17,5022.99,5023.81,5024.33,5024.87,5025.14,5025.46, 
5025.95,5026.22,5026.5,5024.32,5018.87,4995.16,4983.05,4971.65, 
4953.37,4922.25,4913.93,4876.22.4861.85,4857.19,4876.77,4894.18, 
4918.49,4947.27,4973.6,5002.04,5060.12,5089.61, 

NO  = 14,  MU  = 12,  F = 50. , MGSOP  = 4,  WS  = 5. , 

R = 1000. , DELTAR  = 1000. , RMAX  = 26000. , 

$END 


Output  prints  follow  in  Table  7-4. 

UNCLASSIFIED 


Table  7-3.  (U)  Example  2 


The  input  data  deck  for  Example  2 is  presented  below : 
North  Atlantic  Historical  Profile 


$INPUTS 

METRIC  = 1 

IOCEAN  = 2,  IPROFL  - 16,  ISEASN  = 2, 

ZB  = 6500. , 

ZS  = 20. , ZR  = 50. , 

P = 3500.  , MGSOP  = 2,  WS  = 15. , 

R = 500. , DELTAR  = 500. , RMAX  = 80000. , 
IDL  = 47, 


$END 


THEDA(l)  = -60.  ,-29.  ,-28.  ,-26.  ,-25.  ,-23.  ,-22.  ,-21.  ,-20. 
-16. ,-15. ,-14. ,-13. ,-12. ,-11. ,-10.  ,-8. ,-6.  ,-5. ,-3. ,-l.  ,0 
6.  ,8.  ,10.  ,11.  ,12. ,13., 14. ,15.  ,16. ,17.5,18.5,20. ,21. ,22. , 
26. ,28. ,29. ,60. , 


, ,-18.5,-17.5, 
. ,1. ,3.  ,5. , 

23. ,25.  , 


DL(  1)  = 30. ,30. ,20. ,15. ,13. ,11.5,2*11. ,11.5,13. ,15. ,20. 

14..  10.. 8.. 5.. 2. 7. 1.8.. 8.. 2.0... 2.. 8. 1.8. 2. 7. 7. 5.. 8.. 10 

30. . 20. .15. . 13. .11.5.2*11. ,11.5,13. ,15. ,20. ,30. ,30. , 


30. 

,14 


28. ,20.  , 
,20. ,28. , 


Output  prints  follow  in  Table  7-5. 
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TABLE  7-4 

. (U)  Printed  Output* 

for  Example  t 

ii 

92.9667 

1674.7733 

■ e 

12 

100.0000 

1672.9566 

PROFILE  UNITS  (FT. FT/S  OR  DEG 

F): 

13 

133.3333 

1665.0533 

& 

14 

166.6667 

1661.0166 

VELOCITY 

PROFILE: 

15 

200. 0000 

1657.2167 

16 

266.6667 

1651.1233 

N 

Z 

C OR  T 

17 

400. 0000 

1640.7500 

18 

437.4333 

1637.9767 

1 

.0000 

5022.1700 

19 

656.1667 

1625.4066 

it-. 

«-* 

2 

50.0000 

5022.9900 

20 

833.3333 

1620.6166 

3 

100.0000 

5023.8100 

21 

1093.6000 

1619.0633 

4 

131.2000 

5024.3300 

22 

1640.4000 

1625.5900 

5 

164.0000 

5024. 8700 

23 

2187.2333 

1631.3933 

4 

6 

180.5000 

5025.1400 

24 

2734.0333 

1639.4966 

7 

200.0000 

5025.4600 

25 

3280.8333 

1649.0900 

8 

229.7000 

5025.9500 

26 

3827.6333 

1657.8666 

9 

246.1000 

5026.2200 

27 

4374.4333 

1667.3466 

N 

10 

262.5000 

5026.5000 

28 

5468.0666 

1686.7066 

1 1 

278.9000 

5024.3200 

29 

601 4. 8666 

1696.  5367 

12 

300.0000 

5018.8700 

£ 

13 

400.  0000 

4995. 1600 

FREQUENCY  (HZ):  50.00 

14 

500.0000 

4983.0500 

15 

600.  0000 

4971.6500 

MIND  SPEED 

(KTS) : 5.00 

r* 

16 

800.0000 

4953.3700 

17 

1200. 0000 

4922.2500 

SURFACE  LOSS  TABLE  (DEG.DB): 

18 

1312.3000 

4913.9300 

19 

1968.  5000 

4876.2200 

1 

ANGLE 

LOSS 

fl 

20 

2500.0000 

4861.8500 

-• 

21 

3280.8000 

4857.  1900 

1 

.0000 

.0000 

22 

4921.2000 

4876.7700 

2 

2.0000 

.0765 

* 

s 

23 

6581.7000 

4894.1800 

3 

4.0000 

.1542 

.s 

24 

8202.  1000 

4918.4900 

4 

6.0000 

.2332 

25 

9842.  5000 

4947.2700 

5 

8.0000 

.3133 

26 

11482.9000 

4973.6000 

6 

10.0000 

.4311 

Cv 

27 

13123.3000 

5002.0400 

7 

12.0000 

.6388 

28 

16404.2000 

5060.1200 

8 

1 4.0000 

.8578 

29 

18044.6001 

5089.6100 

9 

16.0000 

1.0872 

n 

10 

18.0000 

1.3259 

£» 

SOURCE  DEPTH:  500.00 

11 

20.0000 

1.5729 

12 

22.0000 

1.8275 

^ v 

RECEIVER  DEPTH:  300.00 

13 

24.0000 

2.0890 

v'  ■ 

14 

26.0000 

2.3571 

~ ■*' 

VELOCITY 

PROFILE  (YD. YD/S); 

15 

28. 0000 

2.6317 

16 

30.0000 

2.9130 

N 

DEPTH 

VELOCITY 

17 

32. 0000 

3.2016 

18 

34.0000 

3.4981 

1 

.0000 

1614.0566 

19 

36.0000 

3. 8035 

2 

16.6667 

1674.3300 

20 

38.0000 

4.1 189 

3 

33.3333 

1674.6033 

21 

40. 0000 

4.4457 

4 

43.7333 

1674.7767 

22 

42.0000 

4.7853 

5 

54.6667 

1674.9566 

23 

44. 0000 

5.1392 

W 

6 

60.  1667 

1675.0467 

24 

46.0000 

5.5091 

JR  ■ 

7 

66. 6667 

1675.  1533 

25 

48. 0000 

5.8969 

8 

76„5667 

1675.3167 

26 

50.0000 

6.3044 

V' 

9 

82.0333 

1o75.  4066 

27 

52.0000 

6.  7340 

ri 

10 

87.5000 

1675.5000 

28 

54.0000 

7. 1862 
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29 

56.0000 

7.6698 

15 

’8.0000 

4.8723 

30 

58.0000 

8.1822 

16 

30.0000 

5.2309 

31 

60.0000 

8.7293 

17 

32,  0000 

5.  5666 

18 

34.0000 

5.8804 

MGS  PROVINCE: 

4 

19 

36. 0000 

6.1730 

20 

38.0000 

6.4452 

BOTTOM  LOSS  TABLE  (DEG. 06): 

21 

40.000C 

6.697.' 

22 

42.0000 

6.931  1 

1 

ANGLE 

LOSS 

23 

4/. 0000 

7.1468 

24 

46.0000 

7.3446 

1 

.0000 

.0000 

25 

48.0000 

7. 5254 

2 

2.0000 

.0000 

26 

50.0000 

7.6898 

3 

40  0000 

.0000 

27 

52.0000 

7.8384 

4 

6.0000 

.0000 

28 

54.0000 

7.9716 

5 

8.0000 

.0000 

29 

56. 0000 

8.0899 

6 

1 0.0000 

.4581 

30 

58.0000 

8.1938 

7 

12.0000 

1.0666 

31 

60. 0000 

8.2857 

8 

14.0000 

1.6430 

9 

16.0000 

2.1884 

SOURCE 

ANGLES  (DEG)  FROM  .00  fo  60.00 

10 

18.0000 

2.7040 

1 1 

20.0000 

3.1909 

NORMAL 

MODES  FROM  1 to  10 

12 

22.0000 

3.6502 

13 

24.0000 

4.0829 

REFERENCE  VEIOCITY  CO  (YD):  16/2.9666 

14 

26.0000 

4.4899 
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PROPAGATION  LOSS  TABLE 

: 

50 

50.000 

96.942 

68.445 

51 

51.000 

107.024 

88.553 

1 

R 

PL 

PLRMS 

52 

52.000 

103.411 

88.635 

(KYO) 

<DB) 

(08) 

53 

53.000 

100.820 

88.664 

54 

54.000 

103.274 

88.900 

1 

1.000 

60.244 

59.983 

55.000 

99.064 

88.868 

2 

2.000 

67.918 

68.355 

j 

56,000 

101.841 

88.  886 

3 

3.000 

72. 1 28 

71.714 

57 

57.000 

113.750 

89.116 

4 

4.000 

76. 883 

76.054 

58 

58.000 

108.248 

88.849 

5 

5.000 

79.808 

80.267 

59 

59.000 

101.784 

88.900 

b 

6.000 

81.210 

81.099 

60 

60.000 

94.658 

88.906 

7 

7.000 

82.445 

82.068 

61 

61.000 

92.123 

88.252 

e 

8.000 

84.438 

83.788 

62 

62.000 

88.976 

87.968 

9 

9.000 

85.512 

84.096 

63 

63.000 

85.948 

87.  158 

10 

10.000 

85.479 

85. 105 

64 

64. 000 

85. 093 

84.931 

11 

1 1.000 

96.441 

85. 937 

65 

65.000 

82.432 

83.575 

12 

12,000 

88.825 

86. 109 

66 

66.000 

85.930 

82  . 87  7 

13 

1 3.000 

88.308 

86.984 

67 

67.000 

88.553 

81.989 

14 

14.000 

85.517 

87.376 

68 

68. 000 

82.980 

82.  1 27 

15 

1 5.000 

92.884 

87.492 

69 

69.000 

81.356 

83.228 

1 

16.000 

86.688 

87.959 

70 

70.000 

84.375 

84.430 

17 

17.000 

84.526 

88.066 

71 

71.000 

93.858 

86.501 

18 

18.000 

86.276 

88. 181 

72 

72.000 

97.242 

90.499 

19 

19.000 

1 17.348 

88.388 

73 

73. COO 

93.006 

91.785 

20 

20.000 

93.761 

88. 272 

74 

74.000 

91.158 

91.089 

21 

21.000 

95.390 

88.250 

75 

75.000 

93.436 

92.781 

22 

22.000 

84.059 

88.  176 

76 

76.000 

91.288 

93.336 

23 

23.000 

82.664 

88.039 

77 

77.000 

91.264 

93.063 

24 

24.000 

84.078 

88.043 

78 

78.000 

91.759 

94.250 

25 

25.000 

93.427 

87.959 

79 

79.000 

90.400 

94.  187 

26 

26.000 

96.709 

87.859 

80 

80.  000 

89. 228 

94.  168 

27 

27.000 

89.084 

87.  821 

81 

81.000 

89.279 

94.916 

28 

28.000 

87.510 

87.701 

82 

82.000 

88,834 

94.624 

29 

29.000 

88.601 

87.605 

83 

83.000 

68.096 

94.918 

30 

30.000 

94.661 

87.557 

84 

84.000 

88.364 

95.300 

31 

31.000 

105.940 

87.510 

85 

85.000 

91.661 

95.004 

32 

32.000 

91.731 

87.  546 

86 

86.000 

93.227 

95.238 

33 

33.000 

94.089 

87.567 

87 

87.000 

91.506 

95,253 

34 

34.000 

92.545 

87.473 

88 

88.  000 

90.  878 

95.038 

35 

35.000 

95.792 

87,610 

89 

89.000 

92.006 

95.  182 

36 

36.000 

92.696 

87. 634 

90 

90. 000 

93.822 

95.  107 

37 

37.000 

87.521 

87.659 

91 

91.000 

93.321 

9*. 934 

38 

38. OCO 

84.700 

87. 707 

92 

92.000 

93.128 

95.015 

39 

39.000 

83.219 

87.692 

93 

93.000 

99. 177 

*4.860 

40 

40.000 

82.269 

87.746 

94 

94. 000 

97. 104 

94.900 

41 

41.000 

82.  170 

87  . 848 

95 

95.000 

97.082 

94. 937 

42 

'42.000 

82.134 

87.957 

96 

96.  000 

105.515 

94.857 

43 

43.000 

83.262 

88.127 

97 

97.000 

100.312 

94.967 

44 

44.000 

86.610 

88.225 

98 

98.  000 

99. 131 

94.971 

45 

45.000 

84.412 

88.  290 

99 

99.000 

112.632 

95.017 

46 

46.000 

86.917 

86.418 

100 

100.000 

99.114 

95.098 

47 

47.000 

96.689 

88.420 

101 

101.000 

108.137 

95.039 

48 

48.000 

93.484 

88.455 

102 

’02.  000 

103.610 

95.  268 

49 

49.000 

92.893 

88.509 

103 

103.000 

104.647 

95.448 

UNCLASSIFIED 
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1 04 

104.000 

123.648 

95.  583 

158 

1 58.  OOO 

101.657 

98.  537 

105 

105.000 

107.778 

95.754 

159 

159.000 

102.673 

98.728 

106 

106.000 

111.518 

95.927 

160 

160.000 

100.052 

98. 362 

107 

107.000 

102.244 

95.976 

161 

161.000 

100.792 

93.808 

108 

108. 000 

103. 144 

96. 064 

162 

162.000 

99.717 

99. 1 15 

109 

109.000 

100.907 

96.021 

163 

163.000 

98.949 

99.  198 

110 

110.000 

106.647 

95.959 

164 

164.000 

100.245 

99.312 

111 

111.000 

108.354 

95.897 

165 

165.000 

97.  299 

99.572 

112 

112.000 

101.726 

95.706 

166 

166.000 

100.362 

99.665 

113 

113.000 

95.928 

95.558 

167 

167.000 

99.853 

99.879 

114 

114.000 

95.810 

95. 408 

168 

163.  OOO 

99.574 

100.117 

115 

115.000 

93.756 

95.087 

169 

169.000 

90.943 

100.142 

116 

116.000 

93.578 

94.943 

170 

17C.OOO 

97.458 

100.068 

1 17 

1 17.000 

93.941 

94.708 

1 71 

171.000 

99.  178 

100.078 

118 

118.000 

93.303 

94.4  2 5 

172 

172,000 

99.675 

99.967 

1 19 

1 19.000 

93.047 

94.324 

173 

173.000 

98.800 

99.971 

120 

120.000 

91.788 

94.051 

\V, 

174.000 

103. 176 

99.  744 

121 

121.000 

91.137 

93.919 

175 

(75.000 

100.  160 

99.438 

122 

122.  OOO 

91.765 

93,913 

176 

176.000 

100.464 

99.235 

123 

123.000 

91.403 

93.650 

177 

177.000 

99.759 

98.851 

124 

124.  OOO 

92.538 

93.731 

178 

178.000 

99.617 

98.546 

125 

125.000 

95. 174 

93.687 

179 

179.000 

104.098 

98.202 

126 

1 26.  OOO 

103. 065 

93.209 

180 

180.000 

104.265 

97.834 

127 

127.000 

96.566 

93.249 

181 

181.000 

98.551 

97.612 

128 

129.000 

92.322 

92.776 

182 

192.000 

98. 863 

97.313 

129 

129.000 

93.084 

91.306 

1813 

193.000 

108.471 

97.039 

130 

130.000 

103.267 

90.483 

184 

184.000 

101.533 

96.924 

131 

131.000 

88.797 

89.035 

185 

185.000 

97.310 

96.689 

132 

132.  OOO 

86.368 

87. 550 

186 

186.000 

101.650 

96.  553 

133 

133.000 

91.614 

86.765 

187 

187.000 

108.  736 

96.502 

134 

134.000 

88.432 

86.050 

188 

188.000 

97.378 

96. 327 

135 

135.000 

91.954 

85. 632 

189 

189.000 

99.256 

96.376 

136 

136.  OOO 

83.578 

86.317 

190 

190.000 

106. 502 

96. 299 

137 

137.000 

86.497 

87.398 

191 

191.000 

97.499 

96.  177 

138 

138.  OOO 

87.310 

88.  565 

192 

192.000 

98.074 

96.359 

139 

139.000 

96.461 

91.738 

193 

193.000 

99.699 

95.913 

140 

140.000 

100.355 

95.  564 

194 

194.000 

94.793 

95.587 

141 

141.000 

93.382 

94.315 

195 

195.000 

95,505 

95.337 

142 

142.000 

99.106 

95.269 

196 

196.000 

93.211 

93.735 

143 

143.000 

103.951 

97.950 

197 

197.000 

93.134 

92. 394 

144 

144.  OOO 

96.235 

96.814 

198 

198. OOO 

93.514 

91.341 

145 

145.000 

113.648 

97.519 

199 

199.000 

93.  924 

90.041 

146 

146.  OOO 

104.349 

98.913 

200 

200.000 

94.823 

89.341 

147 

147.000 

103.142 

98.008 

201 

201.000 

95.  180 

88.  ''5? 

148 

148.000 

108.664 

98.652 

202 

202.  OOO 

91.223 

88. 436 

149 

149.000 

98.808 

99.029 

203 

203.000 

88.146 

88.655 

150 

150.000 

103. 033 

98. 305 

204 

204.000 

87.268 

89.656 

151 

151.000 

99.205 

98.  833 

205 

205.000 

89.273 

90.500 

152 

152.  OOO 

101.  101 

98.  664 

206 

206.  OOO 

99.217 

92.  159 

153 

153.000 

102.071 

98.344 

207 

207.000 

99.769 

95.979 

154 

154.  OOO 

104.893 

98.743 

208 

208.  OOO 

95.420 

96. 790 

155 

155.000 

105.136 

98.513 

209 

209.000 

101. 195 

96.147 

156 

156.  OOO 

104.753 

98.494 

210 

210.000 

103.467 

99.  556 

157 

157.000 

104.987 

98.  729 

211 

21 1.000 

101.159 

100.  934 
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212 

212.000 

103.081 

99.228 

213 

213.000 

t06.463 

101.823 

214 

214.000 

103.160 

101.983 

215 

215.000 

108.043 

100.948 

216 

216.000 

104.864 

102.705 

217 

217.000 

107.892 

102.036 

21B 

218.000 

105.668 

101.667 

219 

219.000 

106. 843 

102.763 

220 

220.000 

104.524 

101.951 

221 

221.000 

108.301 

102.184 

222 

222.000 

105.713 

102.886 

223 

223.000 

107.090 

102.369 

224 

224.000 

106.539 

102.921 

225 

225.000 

106.459 

103.221 

226 

226.000 

109.314 

103.086 

227 

227.000 

107.544 

103.879 

228 

226.000 

111.738 

103.968 

229 

229.000 

109.340 

104.174 

230 

230.000 

113.618 

104.925 

231 

231.000 

110.707 

104.560 

232 

232.000 

111.330 

104.897 

233 

233.000 

113.744 

105.416 

234 

234.000 

112.029 

104. 730 

233 

233.000 

114.781 

103.003 

23b 

236.000 

108.601 

104.921 

237 

237.000 

110.377 

104.329 

238 

238.000 

110.  052 

103.806 

239 

239.000 

110.382 

103.354 

240 

240.000 

107.735 

102.331 

241 

241.000 

104.278 

102.114 

242 

242.000 

103.097 

101.435 

243 

243.000 

107.683 

101.057 

244 

244.000 

106.303 

100.743 

245 

245.000 

106.259 

100.278 

246 

246.000 

105,765 

99.946 

247 

247.000 

105.657 

99.581 

248 

248.000 

105. 575 

99.305 

249 

249.000 

103.961 

99.11 1 

230 

250.000 

107. 746 

98.886 

251 

251.000 

106.501 

98.794 

252 

252.000 

106.988 

98. 749 

253 

253.000 

105.224 

98.623 

254 

254. 000 

106. 594 

98.706 

255 

255.000 

107.953 

98.762 

256 

256. 000 

106, 379 

98.818 

257 

257.000 

108.116 

99.062 

258 

258. 000 

108. 846 

99.008 

259 

239.000 

115.082 

99.165 

260 

260. 000 

109.403 

99.144 

TABLE  7-5.  (U)  PRINTED  OUTPUTS  FOR  EXAMPLE  Z 

PROFILE  FROM  HISTORICAL  VELOCITY 
PROFILE  TAPE  FOR  OCEAN  2, 

PROFILE  16.  SEASON  2 


PROFILE  UNITS  <M,  M/S  OR  DEG  C): 


VELOCITY  PROFI  LE: 


N 

2 

C OR  T 

1 

.0000 

1530.0990 

2 

30.4801 

1528.8798 

3 

106.6802 

1523.3934 

4 

152.4003 

1522.7836 

5 

304.8006 

1523.3934 

6 

381.0007 

1524.0030 

7 

457.2009 

1523.3934 

8 

533.4010 

1522.1742 

9 

609.6012 

1519.4310 

10 

685.8013 

1514.8590 

11 

762.0015 

1509.6774 

12 

838.2016 

1504.  8006 

13 

914.4018 

1500.2285 

14 

990.6019 

1 496.  8757 

15 

1066.8021 

1494.1325 

16 

1219.2024 

1491.0845 

17 

1371.6027 

1490.7797 

18 

1524.0030 

1492.3037 

19 

1628.8036 

1495.9613 

20 

2133.  6042 

1499.9227 

21 

2743.2054 

1508.4582 

22 

3657.6072 

1522.4790 

23 

4572.0090 

1537.4142 

24 

5486.4108 

1522.6543 

25 

6400.  81  26 

1568.5039 

SOURCE  DEPTH:  20.00 


RECEIVER  DEPTH:  50.00 


BOTTOM  06 PTH:  6500.00 
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fELCC  ITY 

PROFI 

li£  ( YO,  YD/S): 

15 

28.0000 

8. 1642 

16 

30. 0000 

8.4495 

N 

DEPTH 

VELOCITY 

17 

32.0000 

8.7469 

18 

34. 0000 

9.0564 

1 

.0000 

1673  . 3333 

19 

36.0000 

9 . 3786 

2 

21.8722 

1672.4584 

20 

38. 0000 

9.7138 

3 

33.3333 

1671.9999 

21 

40.0000 

'0.0625 

4 

54.6806 

1670.4629 

22 

42. 0000 

10.4256 

5 

116.6667 

1665. 9999 

23 

44.0000 

10.8039 

6 

166.6667 

1665.3333 

24 

46. 0000 

11. 1984 

7 

333.3333 

1665.9999 

25 

48.0000 

1 1.6106 

8 

416.6666 

1666.6666 

26 

50. 0000 

12.0419 

9 

500. 0000 

1665. 9999 

27 

52.0000 

12.4942 

10 

583.3333 

1664.6666 

28 

54.0000 

12.9699 

11 

666.  6666 

1661.6666 

29 

56.0000 

13.4715 

12 

750.0000 

1656.6666 

30 

58. 0000 

14.0024 

13 

833.3333 

1650.9999 

31 

60.0000 

14.5664 

14 

916.6666 

1645.6666 

15 

1000.0000 

1640.6666 

MGS  PROVINCE: 

2 

16 

1083.3333 

1636.9999 

17 

1166.6666 

1633.9999 

BOTTOM  LOSS  TABLE  (DEG,  DB): 

18 

1333.3333 

1630.6666 

19 

1499.9999 

1630.3333 

1 

ANGLE 

LOSS 

20 

1666.6666 

1631,9999 

21 

1999.9999 

1635.9999 

1 

.0000 

3.5482 

22 

2333.3332 

1640.3333 

2 

2.0000 

3.  8302 

23 

2999.9998 

1649.6666 

3 

4.0000 

4.0597 

24 

3999.9998 

1664.9999 

4 

6. 0000 

4.  2748 

25 

4999.9998 

1681.3333 

5 

8.0000 

4.4773 

26 

5999.9997 

1697.9999 

6 

10.0000 

4.  6684 

27 

6999.9997 

1715.  3333 

7 

1 2.0000 

4.8495 

28 

7108.4722 

1717.2134 

8 

14.0000 

5.  0215 

9 

16.0000 

5. 1853 

FREQUENCY  (HZ):  3500.00 

10 

1 8.  0000 

5. 3416 

1 1 

20.0000 

5.491  1 

WIND  SPEED  (KTS):  15.00 

12 

22. 0000 

5.  6343 

13  v 

24.0000 

5.7718 

SURFACE 

LOSS 

TABLE  (DEG,  DB): 

14 

26.0000 

5.9041 

15 

28.0000 

6.0314 

1 

ANGLE 

LOSS 

16 

30. 0000 

6.  1 54  1 

' 17 

32.0000 

6.2727 

1 

.0000 

6.0013 

18 

34.0000 

6.  3872 

2 

2.0000 

6.0777 

19 

36.0000 

6.4981 

3 

4.0000 

6.  1555 

20 

38.0000 

C.?C.>5 

4 

6.0000 

6.2344 

21 

40.0000 

6.7096 

5 

8.0000 

6.3145 

22 

42.0000 

6.8107 

6 

10. 0000 

6.  3967 

23 

44.0000 

6.9088 

7 

12.0000 

6.4780 

24 

46. 0000 

7.  0043 

8 

14.0000 

6.  581  2 

25 

48.0000 

7.0971 

9 

16.0000 

6.6950 

26 

50. 0000 

7.  1875 

10 

18.0000 

6.9115 

27 

52.0000 

7.2756 

11 

20.0000 

7. 1 392 

28 

54.0000 

7.  3615 

12 

22.0000 

7.3782 

29 

56.0000 

7.4453 

13 

24.0000 

7.6286 

30 

58. 0000 

7. 5271 

14 

26.0000 

7.8906 

31 

60.0000 

7.6070 
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SOURCE  DEVIATION  PATTERN  (DEG,  06)i 


SOURCE  ANGLES  (DEG)  FROM  .00  TO  60.00 


■> 


1 

ANGLE 

LOSS 

NORMAL  MOOES  FROM  1 TO  10 

V 

■ 

1 

-60.0000 

30.0000 

r 

2 

-29.0000 

30.0000 

REFERENCE  VELOCITY  CO  (YD): 

1672.4684 

3 

-28.0000 

20.0000 

4 

-26.0000 

15.0000 

5 

-25. 0000 

13.0000 

v- 

6 

-23.0000 

1 1.5000 

jt 

u 

7 

-22.0000 

11.0000 

8 

-21  .0000 

1 1 .0000 

> 

9 

-20.0000 

11.0000 

10 

-18.5000 

13.0000 

■H 

11 

-17.5000 

15.0000 

*.  M 

12 

-16.0000 

20.0000 

■ 

13 

-15.0000 

30.0000 

s 

14 

-14.0000 

28.0000 

15 

-13.0000 

20.0000 

16 

-12.0000 

14.0000 

•i; 

17 

-11.0000 

10.0000 

18 

-10.0000 

8.0000 

■ 

19 

-8.0000 

5.0000 

V 

20 

-6.0000 

2.7000 

$ 

21 

-5.0000 

1.8000 

22 

23 

-3.0000 

-1.0000 

.8000 

.2000 

24 

.0000 

.0000 

H 

25 

1.0000 

.2000 

s 

26 

3.0000 

.8000 

.V 

27 

5.0000 

1.8000 

28 

6.0000 

2,7000 

l 

29 

8.0000 

5.0000 

30 

1 0.0000 

8.0000 

■ 

31 

11.0000 

10.0000 

32 

1 2.0000 

14.0000 

2 

33 

13.0000 

20.0000 

$ 

34 

14.0000 

28.0000 

3 

35 

15.0000 

30. 0000 

Jv 

36 

16.0000 

20.0000 

E 

37 

17.5000 

15.0000 

l 

38 

18.5000 

1 3.0000 

39 

20.0000 

11.5000 

40 

21.0000 

1 1 .0000 

41 

22.0000 

11.0000 

42 

23.0000 

11.5000 

u 

■ 

43 

25.0000 

13.0000 

44 

26.0000 

1 5.0000 

ft 

45 

28. 0000 

20.0000 

>1 

46 

29.0000 

30.0000 

ft 

47 

60.0000 

30.0000 

s 

m 
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propagation  loss 

TABLE: 

30 

25.000 

11 1.844 

114.901 

31 

25.  500 

122.873 

116.442 

1 

R 

PL 

PLBMS 

52 

26.000 

109.696 

115.804 

(KM) 

(oe> 

(OB) 

53 

26. 500 

118. 599 

114.617 

54 

27.000 

114.568 

115.392 

1 

.500 

32.841 

55.229 

33 

27.  500 

118.749 

113.938 

2 

1.000 

63.009 

62.752 

56 

28.000 

114.270 

114.338 

3 

1.500 

73.006 

68.580 

37 

28.  500 

108.959 

113.492 

4 

2.000 

74.887 

71.01  1 

58 

29.000 

125.304 

113.536 

5 

2.500 

70.126 

68.653 

59 

29.500 

113.781 

113. 882 

5 

3.000 

72.794 

74.590 

60 

30.000 

119.102 

113.225 

7 

3.500 

85.434 

83.255 

61 

30.500 

110.758 

1 1 3.827 

8 

4.000 

85.004 

86.094 

62 

31.000 

122.775 

113.299 

9 

4.500 

87.827 

87.283 

63 

31.300 

11 1.781 

113.586 

10 

5.000 

101.888 

96.066 

64 

32.000 

115.048 

113.633 

11 

5.300 

91.708 

90.423 

65 

32.500 

112.214 

113.517 

12 

6.000 

101.567 

98.332 

66 

33.000 

114.589 

114.159 

13 

6.300 

97.664 

. 94.918 

67 

33.500 

114.984 

113.521 

14 

7.000 

97.191 

96.102 

68 

34.000 

114.244 

114.318 

15 

7.300 

102.036 

100.921 

69 

34.500 

121.091 

114.554 

16 

8.000 

100.142 

97.272 

70 

35.000 

111.959 

114.356 

17 

8.500 

110.164 

104.085 

71 

33.300 

112.534 

114.934 

18 

9.000 

99.730 

99.699 

72 

36.000 

118.674 

114.608 

19 

9.500 

104.843 

104.541 

73 

36.500 

113.298 

115.722 

20 

10.000 

104.143 

106.435 

74 

37.000 

117.770 

115.736 

21 

10.500 

102.462 

102.317 

75 

37.500 

1 19.044 

115.742 

22 

11.000 

110.321 

108.442 

76 

38.000 

118.961 

117.344 

23 

11.500 

101.860 

103.137 

77 

38.500 

117.816 

116.529 

24 

12.000 

108.793 

107.383 

78 

39.000 

121.531 

117.820 

23 

12.500 

105.945 

107.  143 

79 

39.500 

119.142 

118.066 

26 

13.000 

105.232 

105.279 

80 

40.000 

122.882 

118.230 

27 

13.500 

110.371 

112.430 

81 

40.500 

123.228 

1 19.603 

28 

14.000 

105.686 

106.782 

82 

41.000 

118.004 

113.605 

29 

14.500 

112.938 

112.286 

83 

41.500 

120.472 

1 19.700 

30 

15.000 

110.250 

110.829 

84 

42.000 

118.533 

119.179 

31 

15.500 

111.229 

109.823 

85 

42.500 

134.545 

120.333 

32 

16.000 

116.011 

1 16.760 

86 

43.000 

119.870 

121.001 

33 

16.500 

111.723 

109.683 

87 

43.500 

125.900 

119.936 

34 

17.000 

114.831 

116.109 

88 

44.000 

121.789 

121.581 

35 

17.500 

113.620 

112.793 

89 

44.500 

124.359 

120.088 

36 

18.000 

111.719 

113.594 

90 

45.000 

120.244 

121.298 

37 

18.  500 

120.169 

117.009 

91 

45.500 

138.441 

123.451 

38 

19.000 

114.348 

112,407 

92 

46.000 

128.852 

120.240 

39 

19.500 

123.613 

118.251 

93 

46.300 

135.968 

'"■♦.608 

40 

20.000 

116.536 

113.614 

94 

47.000 

134.532 

120.139 

41 

20.500 

114. 146 

113.350 

93 

47.500 

129.309 

121.231 

42 

21.000 

117.137 

1 17.145 

96 

48.000 

123.314 

124.850 

43 

21.300 

114.332 

114.362 

97 

48.500 

126.568 

120.130 

44 

22.000 

126.468 

122.144 

98 

49. 000 

126.312 

125.968 

45 

22.500 

112.762 

115.419 

99 

49.500 

130.204 

1 19.566 

46 

23.000 

121.413 

119.807 

100 

50.000 

132.021 

123.110 

47 

23.500 

116.285 

119.947 

101 

50.500 

130.861 

123.806 

48 

24.000 

1 1 7.  828 

115.483 

102 

51.000 

135.111 

122.119 

49 

24.500 

120.603 

117.924 

103 

51.500 

1 26.667 

125.113 
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104 

52.000 

134.017 

121.011 

109 

52.900 

132.663 

124.651 

106 

53.000 

126.344 

120.131 

107 

53.500 

131.504 

124.102 

106 

54.000 

122.092 

123.054 

109 

54.500 

120.132 

120.645 

110 

55.000 

122.447 

123.095 

111 

55.500 

1 18.976 

119.723 

112 

56.000 

127.042 

123.062 

113 

56.500 

122.117 

124.767 

114 

57.000 

130.615 

1 19.450 

119 

57.500 

121.166 

122.193 

116 

58.000 

121.154 

119.160 

117 

58.  500 

129.443 

119.362 

its 

59.000 

122.415 

121.924 

119 

59.  500 

133.837 

117.097 

120 

60.000 

119.731 

117.819 

121 

60. 500 

124.396 

116. 742 

122 

61.000 

122.668 

1 16.958 

123 

61.500 

121.266 

118.509 

124 

62.000 

123.509 

114.689 

125 

62. 500 

124.831 

119.424 

126 

63.000 

129.829 

1 14. 189 

127 

63. 500 

131.429 

115.048 

128 

64.000 

112.539 

114.331 

129 

64. 500 

97.271 

97. 142 

130 

65.000 

107.940 

97.674 

131 

65.  500 

104.047 

97.757 

132 

66.000 

99.498 

99.245 

133 

66.  500 

101.608 

96.485 

134 

67.000 

103.440 

97.128 

135 

67.  500 

108.530 

100.878 

136 

68.000 

112.609 

100.455 

137 

68.  500 

102.900 

97.650 

138 

69.000 

102.400 

104.268 

139 

69.  500 

103.843 

105.803 
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Appendices  A and  B of  Yarger  (1976)  has 
been  undertaken  to  obtain  the  meaning  of 
each  specific  RAYMODE  product.  The  In- 
clusion of  the  printed  outputs,  there- 
fore, Is  Intended  to  minimise  the  con- 
version process  required  to  adapt  the 
model  at  other  Installations.  References 
to  the  IGS  plotting  system  may  be  re- 
moved by  deleting  the  code  specifically 
marked  for  this  purpose  form  the  listing 
of  the  model's  main  element  RAMODX  and 
deleting  the  four  plotting  subroutines. 

1.  (U)  Print  Option:  IPRINT 

(U)  Optional  prints  of  ray  Information 
are  available  for  each  propagation  model 
Index  J referring  to  the  separate  propa- 
gation types  sustained  by  the  selected 
SVP  in  order  of  Increasing  source  ray 
angle,  where  typically 

J ■ 1 refers  to  surface  duct, 

J - 2 refers  to  convergence  sone, 
and  J - 3 refers  to  bottom  bounce. 

(U)  The  printout  shows  the  velocity  in- 
terval from  CHIN  >Cq  to  CMAX  used  for 
the  propagation  type  and  the  number  of 
trapped  normal  modes  of  that  type.  The 
next  information  reveals  phase  changes 
where  PHIl  is  the  phase  change  undergone 
when  propagation  direction  reverses  from 
an  upward  to  a downward  direction  (ray 
apex  or  surface  reflection),  and  PHI2  Is 
the  corresponding  phase  change  from  a 
downward  to  an  upward  direction  (ray 
nadir  or  bottom  reflection).  The  ray 
path  cycle  limits  used  by  the  program 
for  that  J Index  are  also  printed,  k 
table  of  source  and  receiver  angle  In 
degrees  versus  range  In  yards  and  travel 
time  in  seconds  for  one  cycle  (Q-l) 
and  all  propagation  paths  n from  1 to  4 
at  each  entry  k In  the  ray  table  is  out- 
put along  with  cycle  range  and  cycle 
travel  time.  The  four  paths  between 
source  S and  receiver  R may  be  illus- 
trated using  first  BB(Q  - 1)  as  follows: 


Path  1 Path  2 Path  3 Path  4 


Higher  order  paths  (Q>1)  are  identified 
by  generalising  the  above  scheme.  Short 
range  direct  paths  (Q^>)  consist  of  only 
Path  2 with  Path  1 if  ZS>ZR  or  Path  4 if 
ZS<ZR.  To  compute  horlsontal  range  and 
travel  time  at  cycle  Q-0  (direct  path) 
or  Q>1  for  a particular  path  n 

Ranges  “ (path  n range  for  Q-l)^  + 

(Q-l)  * (cycle  range)k 

Tlrne^  - (path  n time  for  Q“l)jt  + 

(Q-l)  * (cycle  tlme)^ 

for  each  entry  k In  table  Ignoring  nega- 
tive ranges  for  Q-0,  which  restricts  the 
terms  to  the  appropriate  direct  paths. 
The  number  of  points  printed  in  the 
table  will  be  equal  to  NL  or  a value 
selected  by  the  program.  Similar  expres- 
sions to  the  above  will  be  printed  after 
the  table  as  a user  aid  to  hand  calcula- 
tions of  range  and/or  travel  time  at 
particular  source  and/or  receiver 
angles.  This  Information  is  plotted  for 
each  cycle  in  the  angle  versus  range 
plots  and  travel  time  plots  discussed 
below.  To  suppress  or  Include  the  ray 
tracing  outputs  the  user  must  input 
IPRINT  as  follows: 

IPRINT<0  to  cancel  ray  information  print 
56  to  include  ray  information 
print  (default  1) 

2.  (U)  Plot  Options:  PLOTCZ,  PtOTOP,  PLOTT, 
PLOTPl  and  Scale  Term:  PLO 

(U)  The  following  plotting  options  are 
available  from  the  RAYMODE  computer  pro- 
gram with  the  NUSC/NL  IGS  system. 

(1)  SVP  plot  and  a plot  of  bottom 
loss,  surface  loss,  and  deviation  loss 
(beam)  data. 

(2)  Source  and/or  receiver  angle 
versus  horizontal  range  from  each  ray 
path  n from  1 to  4 and  cycle  Q included 
in  propagation  loss  computations. 

(3)  Travel  time  versus  horizontal 
range  for  the  same  ray  paths. 
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(4)  Propagation  loss  versus  horizon- 
tal range  for  coherent  or  random  phase 
combination  of  multlpaths. 

(U)  The  coherent  phase  (PL)  curve  is 
drawn  with  a normal  line  (if  not  dele- 
ted) and  the  random  phase  (PLRMS)  curve 
is  drawn  with  a heavier  line  L(if  not 
deleted).  The  key  containing  frequency, 
source  depth,  and  receiver  depth  again 
appears  to  Identify  the  case  inputs. 

(U)  The  various  user  controls  for  the 
plots  follow , 

PLOTCZ  - integer  plot  option  for  veloc- 
ity profile  and  input  loss 
table  plots,  where 

0 cancels  SVP  and  input  loss 
table  plots. 

> 0 plot  SVP  and  input  loss 
tables  (default  1). 

PLOTOP  - integer  plot  option  for  ray 
angle  versus  range  plots 

£ 0 cancels  ray  angle  versus 
range  plots. 

- 1 plot  only  source  angle  vs. 
range. 

- 2 plot  only  receiver  angle  vs. 
range . 

3 plot  both  source  and  receiver 
angle  vs.  range  on  separate 
graphs  (default  3). 

PLOTT  - integer  plot  option  for  travel 
time  versus  range  plot 

< 0 cancels  travel  time  plot. 

> 0 plot  travel  time  vs.  range 
(default  1). 

PLOTPL  - Integer  plot  option  for  random/ 
coherent  propagation  loss  ver- 
sus range 

- 1 plot  only  coherent  phase 
propagation  loss. 

- 2 plot  only  random  phase  propa- 
gation l08S  • 

■ 1,2  plot  both  coherent  and  ran- 
dom propagation  loss  on  same 
graph  (default  3). 

PLO  - smallest  dB  value  of  propagation 
plot  loss  scale,  PLOX).  (default 
40.). 


8.0  (U)  Organization  Responsibility  for 
RAYMODE  X 

(U)  Responsibility  for  the  RAYMODE  X 
model  resides  at  the  Naval  Underwater 
Systems  Center,  New  London  Laboratory, 
New  London,  Connecticut  06320. 

(U)  For  theory  upon  which  the  model  is 
based: 

Dr.  Gustave  A.  Lelbiger,  Code  327, 
447-4221. 

Mr.  Roy  Deavenport,  Code  321, 
447-4779. 

(U)  For  model  development: 


Dr.  Gustave  A.  Lelbiger, 

Code 

327, 

447-4221. 

Ms.  Donna  F.  Yarger, 

Code 

327, 

447-5198. 

Mr.  Eugene  M.  Smith, 

Code 

327, 

447-5347. 

(U)  For  computer  implementation: 


Ms . Donna  f 

r.  Yarger,  Code 

327, 

447-5198, 

for  UNIVAC. 

Mr . Eugene 

M.  Smith,  Code 

327, 

447-5347, 

for  UNIVAC. 

Mr . * ‘ omas 

Cannan,  Code 

3351, 

447-4738, 

for  Tektronix. 

Mr.  Georpe 

Brown,  Code 

3351, 

447-4680, 

for  PDP  11/70. 

Mr . James 

Bairstow,  Code 

3351, 

447-5514, 

for  HP9845. 

For  model  maintenance:  all  individuals 

listed  above. 

9.0  (U)  Test  Cases  for  implementation  on  a New 
Computer 

(U)  The  only  test  cases  which  can  be 
used  to  check  out  RAYMODE  X on  a new 
computer  are  given  in  Appendices  A and  B 
of  the  User's  Guide  (Yarger,  1976)  and 
also  in  Tables  7-2  and  7-3  (inputs)  and 
Tables  7-4  and  7-5  (printed  outputs). 
The  graphical  outputs  would  be  of  inter- 
est to  only  those  users  with  an  IGS  sys- 
tem. The  second  test  case  utilizes  a 
historical  velocity  profile  (from  a 
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tape)  which  is  also  given  in  the  output. 
A word  of  caution  in  using  these  test 
cases:  the  test  cases  should  be  used 

only  for  the  version  of  RAYMODE  given  in 
Appendix  C of  the  user's  guide.  Many 
versions  of  RAYMODE  exist  which,  due  to 
changes  in  physics,  program  language, 
and  method  of  calculation,  will  not  give 
answers  to  the  test  problems  identical 
to  those  of  the  users's  guide. 

9.1  (U)  Computer  System  on  Which  RAYMODE 
Versions  on  Running 

(U)  Versions  of  RAYMODE  exist  on  the 
following  computers  at  NUSC/NL: 

UNI VAC  1108 

PDP  11/70 

HP  9825 

HP  9845 

Tektronix  4051 

(U)  Outside  NUSC/NL,  many  decks  of  RAY- 
MODE have  been  supplied  to  Navy  contrac- 
tors and  military  laboratories  upon  re- 
quest, as  listed  in  Table  9-1.  Upon  what 
machines  these  decks  have  been  implemen- 
ted Is  not  readily  available.  A poten- 
tial new  user  of  RAYMODE  may,  however, 
recognise  one  of  the  organizations  of 
Table  9-1  as  processing  an  identical 
computer  to  the  one  of  Intended  RAYMODE 
implementation  and,  hence,  gain  from  the 
experience  of  that  organisation.  RAYMODE 
is  known  to  exist  on  IBM  computers  and 
on  the  UYK-20. 

(U)  The  RAYMODE  Univac,  IBM,  and  PDP 
11/70  models  are  written  in  FORTRAN  V; 
the  UYK-20  version  is  in  CMS-2;  the  HP 
9845  and  Tektronix  4051  versions  are  in 
BASIC . 

(U)  The  UNIVAC  versions  of  the  RAYMODE 
model  use  single  precision  arithmetic  on 
a 36-bit  machine.  In  converting  to  a 
computer  system  of  less  precision  (i.e., 
a 32-bit  IBM  360),  revised  versions  of 
subroutine  THREEH,  which  computes  a sen- 
sitive range  derivative,  and  subroutine 
SQUD,  which  computes  a term  based  on  a 
large  machine  dependent  Integer  to 
reduce  the  arguments  of  trigonometric 


functions,  may  be  necessary.  All  UNIVAC 
versions  use  NAMELIST  inputs  which  are 
not  available  in  PDP,  HP,  or  Tektronix 
machines;  these  machines  use  a conversa- 
tional mode  of  input  and  output  parame- 
ter selection  keyed  in  at  the  terminal 
by  the  user,  guided  by  the  software. 

10.0  (U)  RAYMODE  Versions 

(U)  The  following  information  was  sup- 
plied by  NUSC/NL  on  8 August  1980. 

(U)  On  the  UNIVAC  1108  at  NUSC/NL  there 
are  several  current  versions  of  RAYMODE, 
all  in  FORTRAN: 

• The  RAYMODF.  Standard  model,  as  de- 
scribed in  enclosure  (1);  all  previous 
versions  of  the  model  before  1976  (i.e., 
RAYMODE  IV  or  RAYMODE  IX)  are  no  longer 
in  use. 

• The  SUBshA  RAYMODE  version,  as  de- 
scribed by  Weltzner  and  Pearson  (1980); 
this  version  allows  convenient  multiple 
Inputs  (up  to  125)  of  frequency  for 
broadband  computations  and  up  to  ten 
source/receiver  depth  combinations  for  a 
given  environment  and  writes  random 
phase  propagation  loss  vs.  range  outputs 
to  magnetic  tape  or  disc  files. 

e The  RAYMODE  Upgrade  version  being 
implemented  at  System  Consultants,  Inc., 
Arlington,  Virginia,  currently  being 
tested  and  debugged,  contains  a more 
sophisticated  surface  duct  calculation 
and  a revised  "RAYMODE  method"  or  an  FFT 
calculation  as  well  as  more  efficient 
mode  summation  and  bottom  bounce  code; 
this  model  shows  promise  of  significant- 
ly reducing  RAYMODE  run  time  for  some 
cases  while  improving  accuracy. 

• A multi-path  version  of  RAYMODE 
prints  a large  matrix  of  source  angle, 
receiver  angle,  travel  time,  and  propa- 
gation loss  for  each  discrete  range, 
path  cycle,  and  angular  grouping  and  op- 
tionally writes  this  data  to  magnetic 
tape  for  further  processing;  this  model 
is  not  checked  out  at  this  time. 
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Table  9-1.  (U)  List  of  activities  using  RAYMODE  from  1976  to  June  1980 


N Date 


individual 


Organization 


Address 


1 Feb  9,  1976  Ed  Chaika 

2 Jun  10,  1976  David  Michel 

3 Jun  24,  1976  Jerry  Lobdill 

4 Jul  12,  1976  Dr.  Grant  Gartrell 


NUC  (POSSM  Comm.)  San  Diego,  CA  92132 


5 Jul  1976 


Steven  Schuster 


6 Jun  1,  1977  Dr.  Richard  Johnson 

7 Jul  1,  1977  John  Salsbury 

8 Oct  14,  1977  Raymond  R,  Guenther 

9 Dec  , 1978  Bert  Loomis 

LO  Jan  , 1979  Jerry  Bardin 

LI  Mar  21,  1979  Raya  Stern 

L2  Mar  29,  1979  R.J.  Urick 


Magnavox  Co. 


Tracor 

Weapons  Research 
Estab. 

Rockwell  Inter- 
national 

Oregon  State  Univ, 
School  of  Oceano- 
graphy 

NUSC/NPT  Code  444 

Automation  Indus. 
VITRO 

NAVOCEANO 

RADIAN 

Bolt,  Berenak,  & 
Newman 

TRACOR,  Inc. 


L3  May  2,  1979  Dr.  Charles  E.  Schmidt  Honeywell,  Inc. 

Marine  Systems 

L4  May  10,  1979  Bob  Zeskind  BBN 


x/00  Magnavox  Way 
Dept.  529,  Plant  3 
Ft.  Wayne,  IN 

Austin,  TX 


Box  2151  GPO,  Adelaide, 
5A  5001  Australia 

Anaheim,  CA 


Corvallis,  OR  97331 
Newport , RI 

Silver  Springs,  MD  20910 


Austin,  TX 
Cambridge,  MA 

1601  Research  Blvd, 
Rockville,  MD  20850 

5303  Shilshole  Ave,  NW 
Seattle,  Wash.  98107 

1701  Fort  Myer  Drive 
Arlington,  VA  22209 


L5  May  15,  1979  via  Dr.  VonWinkle 


L6  Jul  12,  1979  Hy  Grcenbaum 


L7  Jun  1979 


Ron  Maeur 


L8  Aug  15,  1979  Mike  Turner 


L9  Aug  28,  1979  Ed  Bur ledge 


Robert  McGirr 


SACLANT  ASW 
Research  Center 


Rockwell  Inter- 
national 


TRACOR,  Inc. 

Naval  Ocean 
Systems  Center 


Anaheim,  CA 

7600  Colshire  Drive 
McLean,  VA  22102 

6500  Tracor  Lane 
Austin,  TX  78721 

San  Diego,  CA  92153 
(Code  724) 
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Table  9-1  (cont.).  (U)  List  of  activities  using  RAYMODE  from  1976  to  June  1980 


N 

Date 

Individual 

Organization 

Address 

21 

Oct  19,  1979 

Stewart  Lingley 

System  Planning 
Corp. 

1500  Wilson  Blvd 
Arlingotn,  VA  22209 

22 

Nov  29,  1979 

D.  Paquette 

Code  38211,  NUSC 

NUSC/NPT 

23 

Feb  12,  1980 

Michael  Libby 

IBM,  Manassas 

Bldg  400/041, 

9500  Goodwin  Drive 
Manassas,  VA  22110 

H 

Feb  12,  1980 

W.H.  Lunceford 

Naval  Training 
Equip.  Center 

Code  N233 
Orlando,  FL  32813 

25 

Mar  5,  1980 

LCdr  Frank  Hiestand 

C0M0PTEVF0R  Naval 
Base 

Norfolk,  VA  23511 

Mar  1980 

ASA  Davis 

NUSC 

NUSC/NPT 

27 

Apr  5,  1980 

George  E.  Miller 

RCA  Automated 
Systems 

PO  Box  588  MS  1-1 

Burlington,  MA  01803 

28 

Mar  28,  1980 

Joe  Fenier 

McDonald  Dougas 
Aeronautics 

PO  Box  516 

St.  Louis,  MO  63166 

29 

Apr  2,  1980 

Bob  Woodham 

Naval  Surface 
Weapons  Center 

Code  031 

Silver  Spring,  MD  20910 

May  10,  1980 

Paul  Scherer 

Naval  Air  Develop- 
ment Center 

Warminster,  PA  18974 

31 

Jun  19,  1980 

Stephen  P.  Koch 

B-K  Dynamics,  Inc. 

15825  Shady  Grove  Road 
Rockville,  MD  20850 
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(U)  On  other  computers  within  NUSC/NL, 
there  exists: 

s An  HP  9845  version  of  RAYMODE  writ- 
ten In  BASIC  language.  Balrstovr  and 
Medieras  (1980)  describe  this  version 
and  give  a test  case. 

s An  HP  9875  version  also  in  BASIC 
used  by  developer  Dr.  G.  A.  Leibiger  for 
testing  new  concepts  and  techniques; 
this  version  has  inputs  keyed  in  by  user 
and  has  limited  graphical  output. 

e A Tektronix  4051  version  also  in 
BASIC,  again  used  by  Dr-  Leibiger,  as  an 
experimental  model. 

s A PDP  11/70  version  used  by  the  Sur- 
face Ship  Department  written  in  FORTRAN. 

• A Tektronix  4051  version  as  imple- 
mented from  the  propagation  loss  tape  of 
the  Fleet  Mission  Library  using  conver- 
sational input  guidance  and  selectable 
varied  output  graphics  used  in  the 
fleet. 

(U)  A ROM  version  of  RAYMODE  for  the 
Tektronix  is  nearly  ready,  which  is 
estimated  to  improve  RAYMODE  run  time  by 
a factor  of  ten. 

11.0  (U)  Test  Cases  Used  in  the  Evaluation  of 
RAYMODE  X 

(U)  Test  cases  were  chosen  from  experi- 
mental data  sets.  The  experimental  data 
sets  are  described  in  detail  by  Martin 
(1982)  and  constitute  a Portable  Test 
Package  for  model  evaluation.  A subset 
of  these  cases  were  selected  for  the 
RAYMODE  X evaluation  based  on  time  and 
cost  constraints  and  availability  of  the 
data  during  this  evaluation.  The  experi- 
mental sets  selected  were  SUDS,  HAYS- 
MURPHY,  PARKA  II,  BEARING  STAKE,  JOAST, 
LORAD,  FASOR,  and  GULF  OF  ALASKA.  Some 
general  information  on  these  data  sets 
is  given  in  Table  11.1.  As  can  be  read 
from  this  table,  the  data  sets  were 
selected  to  provide  broad  geographic  and 


frequency  coverage  and  coverage  with  re- 
dundancy of  the  various  propagation  mod- 
els. Specific  characteristics  of  the 
subsets  selected  for  the  RAYMODE  X eval- 
uation are  given  in  Table  11.2a-h  which 
include  source  and  receiver  depths, 
mixed  layer  depth  and  depths  of  sound 
channel  axis  and  bottom,  frequency,  and 
maximum  range  of  data.  Sound  speed  pro- 
files, bottom  loss  versus  grazing  angle 
curves,  and  the  measured  acoustic  data 
are  found  in  Appendices  IIA  - IIH,  re- 
spectively, for  the  experimental  sets 
mentioned  above.  The  bottom  loss  versus 
grazing  data  is  that  associated  with  the 
RAYMODE  X model  with  bottom  loss  deter- 
mined from  geographic  area  designator 
charts  with  the  following  exceptions: 

(1)  for  PARKA  II  and  HAYS -MURPHY,  FNOC/ 
N00  bottom  loss  versus  grazing  angle 
curves  used  by  the  FACT  PL9D  model 

(evaluated  by  AMEC  as  reported  in  Volume 
II  of  this  series)  with  bottom  loss 
being  determined  from  NAVOCEANO  area 
designator  charts  was  used  for  RAYMODE  X 
in  addition  to  the  standard  RAYMODE  in- 
puts; (2)  a constant  bottom  loss  of  50 
dB  was  used  for  all  SUDS  cases,  effec- 
tively eliminating  bottom  interacting 

paths,  since  the  experiment  resulted 
from  pulBed  transmission  and  bottom 
reflected  paths  were  time-gated  out;  and 
(3)  bottom  loss  measurements  from  the 
experiment  site  were  used  for  BEARING 
experiment  site  were  used  for  BEARING 
STAKE  since  they  differed  so  greatly 

(i.e.,  show  much  less  loss)  than  would 
be  obtained  from  either  MGS  or  NAVOCEANO 
(used  for  FACT)  area  designator  charts 
and  their  associated  curves. 

11.1  (U)  Results  of  Tost  Coses  Used  in  the  AMEC 
Eviluitwn  of  RAYMODE  X 

(U)  SUDS:  (1)  In  terms  of  decibel  dif- 
ferences between  SUDS  and  RAYMODE  data, 
only  Cases  II,  IV,  and  VI  show  fair 
agreement.  Cases  II  and  IV  are  cases  of 
cross-layer  source/receiver  geometry  at 
200  and  1000  Hz,  respectively.  Cases  of 
cross-layer  geometry  at  higher  frequen- 
cies shows  basic  qualitative  and  quanti- 
tative disagreement  between  SUDS  data 
and  RAYMODE  X results.  All  cases  of  both 


100 


CONFIDENTIAL 


5 


ft 


;$j 

wi 


x 1 


iS 


k.T 


.v 


& 

.• : 

1 00 


V> 

v> 


££! 

& 


TABLE  11. 2A.  (U)  SUDS:  TEST  CASE  CHARACTERISTICS 
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TABLE  11,28.  HAYS-MURPHY:  TEST  CASE  CHARACTERISTICS  (U) 


CASE 

SOURCE 

DEPTH 

(m) 

RECEIVER 

DEPTH 

(m) 

FREQUENCY 

(Hz) 

SOUND 

AXIS 

DEPTH  (m) 

BOTTOM 

DEPTH 

(m) 

I 

24.4 

137.2 

35.0 

61 

2750 

II 

24.4 

137.2 

67.5 

61 

2750 

III 

24.4 

137.2 

100.0 

61 

2750 

IV 

24.4 

137.2 

200.0 

61 

2750 

V 

24.4 

106.7 

35.0 

61 

2750 

VI 

24.4 

106.7 

100.0 

61 

2750 

CONFIDENTIAL 


TABLE  11. 2C.  PARKA:  TEST  CASE  CHARACTERISTICS  (U) 


CASE 

FREQUENCY 
( #2  ) 

SOURCE 
DEPTH 
( m ) 

RECEIVER 
DEPTH 
( m ) 

LAYER 
DEPTH 
( m ) 

SOUND 

AXIS 

DEPTH  ( m } 

BOTTOM 
DFPTH 
( 'm  ) 

I 

50 

152.4 

91.4 

80 

1000 

5500 

II 

400 

152.4 

91.4 

80 

1000 

5500 
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TABLE  11. 2D.  BEARING  STAKE:  TEST  CASE  CHARACTERISTICS  (U) 
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source  and  receiver  In  the  layer  ahov 
lack  of  agreement  between  SUDS  and  RAY- 
MODE  . Of  the  two  below-layer  cases,  VI 
showed  good  agreeaent  between  SUDS  and 
RAYMODE,  but  VIII  did  not.  (2)  The  fig- 
ure of  aerlt  analysis  shows  RAYMODE 
predictions  of  detection  range  to  be 
strongly  pessiaistlc  compared  to  SUDS 
results  in  Cases  I,  II,  III,  IV,  V,  IX 
and  X (not  uncommonly  by  a factor  of 
two) . Only  in  Case  VIII  were  RAYMODE 
detection  ranges  much  longer  (by  a fac- 
tor of  two)  than  those  of  SUDS.  No  clear 
trend  emerges  with  regard  to  source/re- 
celver  geometry.  (3)  It  would  appear 
from  comparison  with  SUDS  data  that  the 
surface  duct  aodule  of  RAYMODE  X is  de- 
ficient. This  is  further  borne  out  in 
Section  3.0,  The  Physics  of  the  RAYMODE 
X Model,  by  R.  Deavenport. 

(U)  HAYS -MURPHY;  (1)  Significant  differ- 
ences in  mean  levels  were  primarily 
responsible  for  pessimistic  detection 
range  predictions  by  the  RAYMODE.  These 
differences  appear  to  be  attributable  to 
the  bottom  loss  inputs  for  the  first 
bottom  bounce  region  (l.e.,  to  25  km). 
Beyond  this  range,  differences  are  as 
great  and  unexplained,  but  bottom  loss 
is  not  a factor.  It  is  to  be  noted  that 
for  this  scenario,  RAYMODE  and  FACT  bot- 
tom loss  inputs  led  to  essentially  the 
same  results.  (2)  haan  differences  be- 
tween the  Hays -Murphy  data  and  the  RAY- 
MODE model  were  smaller  by  about  2 dB 
for  incoherent  results  as  compared  to 
coherent  results.  This  is  reversed  for 
standard  deviations  of  differences  be- 
tween the  model  and  Hays-Murphy  data, 
for  which  RAYMODg  coherent  generally 
showed  about  2 dB  greater  standard  devi- 
ation than  did  RAYMODE  incoherent.  The 
net  effect  is  that  the  RAYMODE  incoher- 
ent curve  is  in  better  agreement  with 
the  experimental  data  than  is  the  RAY- 
MODE coherent  curve  with  regard  tc  gen- 
eral characteristics  (i.e.,  shape),  but 
the  RAYMODE  coherent  results  are  more  in 
agreement  with  the  experimental  data 
with  regard  to  detection  range  coverage 
(although  the  agreement  is  far  from  sat- 
isfactory) . This  is  understandable  since 
detections,  particularly  of  a zonal 


nature,  are  determined  more  by  fluctua- 
tions than  by  mean  levels  (particularly 
for  average  slgnal-to-noise  ratios  which 
are  negative  or  near  zero) . 

(U)  PARKA:  (1)  The  FNOC  bottom  loss 
used  in  RAYMODE  leads  to  better  agree- 
ment with  PARKA  data  than  does  use  of 
the  MGS  bottom  loss.  (2)  At  50  Hz,  the 
three  convergence  zones  as  predicted  by 
RAYMODE  are  successively  more  severely 
retarded  in  range  compared  to  the  PARKA 
data.  The  peak  levels  of  the  RAYMODE 
first  and  second  zones  ere  2 and  1 dB 
less  than  PARKA* a peak  levels,  respec- 
tively. The  third  zone  RAYMODE  peak  is  4 
dB  less  than  PARKA's.  (3)  At  50  Hz,  the 
PARKA  data  shows  substantially  less  loss 
than  does  RAYMODE  in  all  bottom  bounce 
regions  (using  RAYMODE's  MGS  values);  at 
400  Hz  the  trend  is  the  same,  but  agree- 
ment is  somewhat  closer.  (4)  At  50  Hz 
the  PARKA  data  and  RAYMODE  predictions 
agree  well  (near  the  RAYMODE  coherent 
peaks)  in  the  first  bottom  bounce  region 
and  agree  well  in  the  middle  of  the  sec- 
ond bottom  bounce  region.  In  the  third 
bottom  bounce  region,  the  RAYMODE  pre- 
diction shows  less  loss  than  the  PARKA 
data.  This  is  in  contrast  to  the  400  Hz 
results  for  'triiich  PARKA  showed  signifi- 
cantly less  loss  in  all  bottom  bounce 
regions  compared  to  RAYMODE.  (5)  At  50 
Hz  the  first  convergence  zone  as  -e- 
dicted  by  RAYMCDE  is  in  very  good  agree- 
ment with  that  of  PARKA,  but  is  slightly 
wider.  The  second  RAYMODE  CZ  is  found  at 
shorter  range  than  is  PARKA's  by  about 
5 km.  This  situation  is  exaggerated  in 
the  case  of  the  third  convergence  zone. 
At  400  Hz  the  results  for  the  first  and 
second  CZ  starts  are  basically  the  same. 
The  RAYMODE  first  CZ  is  wider  than 
PARKA's,  and  the  second  is  narrower. 

(U)  BEARING  STAKE : (1)  RAYMODE  X coher- 
ent predictions  are  in  better  agreement 
with  BEARING  STAKE  data  than  are  RAYMODE 
X incoherent  predictions.  (2)  Agreement 
between  RAYMODE  and  BEARING  STAKE  re- 
sults are  often  in  reasonable  agreement 
to  ranges  from  as  far  as  60  to  150  km. 
(3)  In  the  difference  curves,  there  is 
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an  underlying  trend  causing  the  differ- 
ence between  BEARING  STAKE  and  RAYMODE 
results  to  become  Increasingly  negative 
with  range.  This  suggests  that  a higher 
critical  angle  In  the  bottom  loss  versus 
gracing  angle  curve  would  lead  to  better 
agreement,  offsetting  this  trend.  (4) 
BEAP'NG  STAKE  data  for  the  receiver  on 
and  *0  m off  the  bottom  show  strong  In- 
terference patterns.  RAYMODE  X predic- 
tions show  patterns  which  are  generally 
out  of  phase  with  those  of  BEARING  STAKE 
at  short  ranges  < <30  km)  and  are  dissim- 
ilar at  longer  ranges.  The  RAYMODE 
interference  patterns  are  generally 
stronger  (l.e.,  greater  peak-to-peak 
excursion)  than  are  those  of  BEARING 
STAKE  data.  (5)  Detection  coverage  re- 
sult 1 are  usually  In  rough  agreement  for 
figures  of  merit  of  75  and  80  dB  between 
RAYMOEE  X predictions  and  BEARING  STAKE 
data.  This  agreement  often  extended  to 
85  and  90  dB  and  in  one  case  to  95  dB. 
(6)  For  FOM  > 95  dB,  BEARING  STAKE  de- 
tection coverage  was  to  much  longer 
range  and  was  more  complete  (l.e.,  bet- 
ter percentage  coverage)  than  were 
RAYMODE's. 

(U)  LORAD:  For  Cases  lA-iG  for  which  the 
source  in  the  surface  duct  and  the  re- 
ceiver at  the  duct  limit  of  100  feet 
(30.5  m)r  (1)  LORAD  shows  somewhat  bet- 
ter detection  coverage  than  RAYMODE  X in 
the  bottom  bounce  regions.  Changing  the 
bottom  class  in  the  model  should  not 
provide  significant  Improvement,  ' since 
the  shapes  of  the  propagation  loss  ver- 
sus range  curves  for  LORAD  and  RAYMODE 
are  different  ir  the  bottom  bounce  re- 
gions. (2)  The  L0R4D  first  convergence 
zone  has  one  lobe  whereas  RAYMODE  X'e 
has  two.  (3)  The  LORAD  C 7.  start  precedes 
RAYMODE's  in  the  first,  third,  and  fifth 
CZs  and  the  converge  is  true  for  the 
second,  fourth,  and  seventh  CZs.  (4;  The 
ends  of  the  LORAD  CZs  occur  at  equal  or 
greater  ranges  than  do  RAYMGDL’s.  (5) 
Fluctuations  in  the  convergence  zones 
are  greater  for  LOR.  r data  than  for  RAY- 
MODE predictions.  (6)  The  RAYMODE  X 
coherent  results  are  in  better  agreement 
with  LORAD  data  than  are  RAYMODE  X inco- 
herent results. 


(U)  For  Cases  IIA-IIG  for  which  the 
source  is  in  the  surface  duct  and  the 
receiver  is  below  the  duct  at  1000  feet 
(305  m):  (1)  LORAD  data  and  RAYMODE  pre- 
dictions differ  by  about  10  dB  in  bottom 
bounce  regions,  with  LORAD  showing  less 
propagation  loss.  The  use  of  lower  bot- 
tom loss  as  model  input  would  lead  to 
better  agreement.  (2)  The  start  and  end 
ranges  for  LORAD  and  RAYMODE  agree  for 
the  first  and  second  convergence  zones 
(averaging  over  all  figures  of  merit). 
The  first  CZ  is  double-lobed  for  both 
model  and  experimental  data.  (3)  From 
the  third  CZ  on,  the  start  of  the  LORAD 
CZ  precedes  that  of  RAYMODE  and  the 
model's  CZ  start  is  steeper  than  that 
indicated  by  the  experimental  data. 
Range  disparities  as  great  as  10  km  are 
found  in  CZ  start  range  between  LORAD 
and  RAYMODE.  (4)  The  range  at  which  the 
convergence  zone  ends  for  LORAD  data  is 
equal  to  or  greater  than  that  for  RAY- 
MODE by  as  much  as  6 km.  (6)  Fluctua- 
tions in  the  LORAD  data  are  of  roughly 
the  same  magnitude  as  the  unsmoothed 
RAYMODE  coherent  output.  (7)  The  RAYMODE 
X coherent  prediction  is  generally  in 
better  agreement  with  the  LORAD  data 
than  is  the  RAYMODE  X incoherent  predic- 
tion. 

(U)  JOAST:  RAYMODE  X predictions  and 
JOAST  experimental  data  were  In  basical- 
ly good  agreement  with  regard  to  conver- 
vergence  zone  shape,  peak  level,  start 
range,  and  zone  duration.  Exceptions  are 
(1)  in  Coses  II  and  III  the  RAYMODE  CZ 
end  Mas  broader  by  about  2-3  km  than 
that  of  JOAST,  and  (2)  the  RAYMODE  and 
JOAST  CZs  were  displaced  In  level  by  10 
dB  (JOAST  exhibiting  greater  loss)  In 
Cases  XII-XIV.  The  use  of  FN0C  bottom 
lots  may  have  been  responsible  for  the 
greater  evidence  of  bottom  loss  contami- 
nation of  the  convergence  zone  (particu- 
larly at  its  end)  for  the  RAYMODE  pre- 
dictions as  compared  to  JOAST  data.  This 
is  true  for  station  3 where  relatively 
less  loss  was  predicted  from  the  FN0C 
bottom  loss  charts  (FN0C  Type  3 for  sta- 
tion 3).  This  compares  to  MGS  Type  6 for 
station  3.  This  is  not  true  for  stations 
1 and  2 where  both  are  characterized  by 
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“J  FNOC  Type  2 and  MGS  Type  2 (Note:  The 

MGS  and  FNOC  bottom  loss  vs.  grazing 
I angle  carves  differ,  but  by  less  than  1 

V dB  over  the  full  angular  extent). 


(U)  FASOR : Aside  from  Case  Ila,  agree- 
ment between  FASOR  data  and  RAYMODE  X 
predictions  Is  qualitatively  good.  Far 
Case  Ila,  the  model  Is  utilizing  a lov 
loss  bottom  and  a significant  reduction 
of  levels  Is  probably  unachievable  with 
available  bottom  loss  curves.  Overall, 
there  Is  a slight  tendency  for  the  FASOR 
results  to  show  less  loss  and  give 
slightly  better  detection  coverage,  al- 
though this  varies  from  case  to  case  and 
also  depends  upon  the  RAYMODE  coherence 
option  chosen. 
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<U)  GULF  OF  ALASKA:  (1)  Fluctuations 

were  not  as  rapid  for  RAYMODE  coherent 
as  for  Gulf  of  Alaska  (GOA)  data  and 
were  larger  for  GOA  by  factors  as  great 
as  two,  with  the  exceptions  of  Case  I 
where  fluctuations  were  10-15  dB  for 
each.  Case  XI  with  3 dB  fluctuations. 
Cases  XII  and  XIII  with  5-15  dB  fluctua- 
tions and  Case  XIV  where  GOA  fluctua- 
tions were  12  dB  and  RAYMODE* a varied 
between  5 and  20  dB.  Note:  the  final 
four  cases  are  in  the  sound  channel  with 
the  source  below  the  axis;  frequency  of 
2.5  kHz.  (2)  The  1.5  kHz  results  (Cases 
I-VI)  show  basic  similarity  between  GOA 
and  RAYMODE  coherent  results.  The  RAY- 
MODE Incoherent  curve  usually  provided  a 
low  propagation  loss  envelope  for  the 
un8moothed  GOA  data.  (3)  There  Is  good 
agreement  In  start  and  end  ranges  be- 
tween GOA  and  RAYMODE  convergence  zones 
except  that  In  Case  VII,  a CZ  was  pre- 
dicted by  RAYMODE  Incoherent  but  not  by 
RAYMODE  coherent,  nor  was  it  evident  in 
GOA  date.  (4)  Results  for  the  1067  m 
source  were  Inconsistent.  More  often 
than  not,  however,  RAYMODE  results  cor- 
responded to  low-loss  envelopes  of  GOA 
data  (Case  VII  showed  severe  unexplained 
disagreement).  (5)  For  the  last  four 
cases  (305  m source)  the  RAYMODE  curves 
were  generally  low-loss  envelopes  for 
GOA  data  or  were  parallel  to  the  low- 
loss  envelope  but  with  3 dB  less  loss. 
(6)  For  the  1.5  kHz  data,  figure  of 
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merit  analysis  shows  the  extent  of  range 
coverage  of  GOA  data  to  be  greater  than 
or  equal  to  that  of  RAYMODE  output.  This 
trend  is  reversed  for  the  2.5  kHz  data 
for  which  RAYMODE  had  consistently  bet- 
ter range  coverage  than  did  GOA  data. 
(7)  Over  range  Increments  where  both 
RAYMODE  and  GOA  data  had  detection  cov- 
erage, the  zonal  detection  coverage 
(detections  per  opportunity  In  percent) 
was  greatest  for  RAYMODE  incoherent 
(usually  100%)  and  least  for  GOA  data 
due  to  fluctuations.  Two  exceptions  were 
Case  II  for  which  the  ZDCs  were  about 
equal  and  Case  IV  for  tftiich  the  GOA  data 
had  better  ZDC  than  either  RAYMODE  op- 
tion. 

12.0  (U)  Summary  and  Racommandations 

(U)  The  RAYMODE  model  produces  propaga- 
tion loss  as  a function  of  range  and 
frequency  In  an  environment  character- 
ized by  a single  sound  speed  profile  and 
a horizontal  ocean  floor.  The  evaluation 
herein  reported  has  been  for  a specific 
version  of  the  RAYMODE  model,  which  Is 
known  as  RAYMODE  X,  and  all  test  cases 
were  :un  on  the  UNIVAC  1108  computer 
(wit«.  EXEC  VIII  compiler)  at  the  New 
London  Laboratory  of  the  Naval  Under- 
water Systems  Center.  Although  RAYMODE 
has  been  primarily  developed  for  tacti- 
cal sonar  applications,  It  has  been  used 
over  the  span  of  frequencies  Implied  by 
surveillance  through  torpedo  explica- 
tions. RAYMODE  can  produce  both  coherent 
and  incoherent  transmission  loss 
results;  however,  when  used  for  fleet 
applications  (as  opposed  to  research 
usage),  the  incoherent  phase  addition 
option  Is  selected.  The  fleet  user  is 
not  given  the  option  of  selecting  values 
for  various  program  controls.  These  are 
minimum  and  maximum  number  of  ray  cycles 
for  all  propgation  modes,  maximum  number 
of  ray  cycles  for  bottom  bounce,  posi- 
tive and  negative  minimum  sonar  angle 
(default  0°),  positive  and  negative  max- 
imum sonar  angle  (default  60°),  the 
first  relative  normal  mode  processed  by 
mode  summation  (default  1),  the  last 
relative  normal  mode  processed  by  mode 
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summation  (default  0),  the  maximum  num- 
ber of  modes  processed  by  mode  summation 
(default  10),  and  the  number  of  points 
in  ray  tables  (default  10).  A detailed 
description  of  these  program  controls  is 
given  in  section  7.0.  The  choice  of 
these  parameters  can  have  a large  effect 
on  running  time.  In  all  cases  used  in 
this  evaluation,  default  values  were 
assigned  for  these  parameters  with  one 
exception  (in  one  FASOR  case,  the  maxi- 
mum sonar  angle  default  of  60°  led  to 
unrealistically  high  values  of  propaga- 
tion loss  over  a short  range  extent; 
changing  the  maximum  sonar  angle  to  85° 
eliminated  the  problem).  Under  default 
conditions,  run  times  on  the  UNIVAC  1108 
varied  from  extremes  of  5.1  to  54.1  sec- 
onds (usually  between  6 and  30  seconds) 
and  were  strongly  scenario  dependent. 
The  number  of  range  points  at  which 
propagation  loss  was  calculated  in  these 
runs  varied  from  200  to  400.  It  should 
be  noted  that  by  altering  the  default 
values  of  the  program  controls,  in- 
creased accuracy  can  be  achieved  at  the 
expense  of  running  time. 

(U)  The  computer  core  required  by  RAY- 
MODE  X (when  dimensioned  for  400  ranges 
and  50  modes  and  ray  points)  without 
plots  in  17319  decimal  words.  When 
dimensioned  for  200  ranges  and  25  modes 
and  ray  points,  the  core  required  (with- 
out plot 8)  is  15894  decimal  words.  We 
recall  that  the  defaults  for  numbers  of 
modes  and  ray  points  is  10. 

(U)  At  NUSC/NLL,  versions  of  RAYMODE 
exist  on  the  following  computers:  UNIVAC 
108,  PDP  11/70,  HP  9825,  HP  9845,  and 
Tektronix  4051.  Of  these,  the  UNIVAC  and 
PDP  versions  are  written  in  FORTRAN  V, 
and  the  HP  and  Tektronix  versions  in 
BASIC.  A version  also  exists  on  a UYK- 
20  computer  written  in  CMS-2  language. 
RAYMODE  versions  support  the  following 
fleet  systems  and  applications:  TRIDENT 
Optimum  Mode  Selection  (0MS),  BQQ-5  so- 
nar (MS,  BQQ-6  sonar  0MS,  Improved  Souar 
Processing  Equipment  (ISPE),  Submarine 
Active  Detection  System  (SADS),  Subma- 
rine System  Effectiveness  and  Assessment 


(SUBSEk j , Sonar  In-Situ  Mode  Assessment 
System  (SIMAS),  and  Fleet  Mission  Li- 
brary. 

(U)  The  physics  of  the  RAYMODE  X model 
was  examined  by  Roy  Deavenport  of  NUSC/ 
NLL,  New  London,  Conn.  The  resultant 
description  (c.f.,  section  3.0)  repre- 
sents the  most  detailed  description  of 
RAYMODE 's  known  physics.  Indeed,  a pri- 
mary deficiency  of  RAYMODE  is  its  lack 
of  documentation  both  in  the  form  of 
reports  describing  the  physics  of  the 
model  and  its  implementation  on  a sub- 
routine-by-subroutine  basis  and  in  the 
form  of  comment  cards  within  the  RAYMODE 
computer  code.  Deavenport  selected  two 
test  cases  for  RAYMODE.  The  first  case 
tested  surface  duct  propagation  with  the 
result  that  RAYMODE  can  yield  poor  sur- 
face duct  results  at  the  lower  frequen- 
cies (<200  Hs).  The  second  case  was  con- 
cerned with  low-frequency  propagation  in 
an  environment  wnere  the  sound  speed 
profile  has  a depressed  channel  above  a 
deep  sound  channel.  It  was  concluded 
that  RAYMODE  does  not  properly  account 
for  depressed  channel  propagation  at  the 
lower  frequencies  (<300  Hz).  This  ap- 
pears to  be  related  to  the  fact  that  in 
RAYMODE  there  is  no  consideration  of 
partial  trapping  of  energy  within  the 
depressed  channel. 

(U)  Two  test  cases  for  RAYMODE  X exist 
in  the  RAYMODE  User's  Guide  (Yarger, 
1971)  but  are  not  necessarily  applicable 
to  other  versions  of  RAYMODE.  In  fact, 
it  appears  that  RAYMODE  is  not  under 
configuration  management.  It  is  acknowl- 
edged that  different  versions  have 
slightly  different  physics  in  some  sub- 
routines (this  Includes  the  surface  duct 
treatment),  are  programmed  in  differe.it 
languages,  and  are  run  on  different  com- 
puters with  different  word  lengths. 
There  do  not  appear  to  be  test  cases  for 
which  all  versions  have  been  run  to  as- 
sure uniformity  of  results. 

(U)  Responsibility  for  the  RAYMODE  X 
model  resides  at  the  Naval  Underwater 
Systems  Center,  New  London  Laboratory, 
New  London,  CT  06320,  for  the  model's 
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theory,  development,  computer  implom^i. 
tatlon  and  model  maintenance  aa  detailed 
In  section  8.0. 

(U)  Results  of  test  cases  (see  sections 
11.0  and  11.1)  Indicate  (a)  basic  dis- 
agreement between  SUDS  data  and  RAYMODE 
X output  for  surface  duct  propagation, 
(b)  basic  agreement  of  convergence  cone 
levels  and  range  for  PARKA  and  JQAST 
data  with  somewhat  less  agreement  with 
LORAD  data,  (c)  variable  agreement  in 
bottom  bounce  regions  with  best  agree- 
ment between  RAYMODE  and  experimental' 
data  found  in  the  first  bottom  bounce 
region  with  successive  deterioration  of 
agreement  as  successively  distant  bottom 
bounce  regions  are  encountered. 

(U)  Useful  features  present  in  RAYMODE  X 
are: 

e user  specification  of  source  and  re- 
ceiver beam  patterns, 

e availability  of  eigenray  information 
Including  travel  times, 

e Independent  specification  of  initial 
range  of  range  Increment,  and 

• user  specification  of  bottom  loss 
versus  gracing  angle  table. 

(U)  The  variety  of  RAYMODE  versions 
available  implies  that  different  results 
are  obtainable  from  all  models  carrying 
the  name  RAYMODE  for  the  same  environ- 
mental inputs  and  the  same  program  con- 
trol selections. 

(U)  In  light  of  the  above,  the  following 
recommendations  are  given: 

• Improve  RAYMODE  documentation,  spe- 
cifically describing  the  physics  and  its 
implementation  on  a subroutine-by-sub- 
routine  basis,  defining  all  physical 
variables  and  giving  their  FORTRAN  (aljo 
Basic  and  CMS-2)  names,  providing  flow 
charts  and  internal  to  the  program,  pro- 
viding adequate  comment  cards  to  allow 
the  user  to  trace  the  program  logic. 


s Initiate  a "fault-finding”  evalua- 
tion of  RAYMODE  physics  whereby  test 
cases  may  be  designed  t valuate  spe- 
cific aspects  of  RAYMODE  ysics  includ- 
ing appro  -imations  and  as  options. 

t Test  the  effect  of  word  length  on 
the  itAYMODE  X output.  Determine  if  and 
where  the  use  of  double  precision  arith- 
metic is  indicated. 

• Provide  a variety  of  test  cases  to 
assure  that  transfer  of  RAYMODE  X to  a 
new  computer  is  accomplished  with  RAY- 
MODE providing  consistent  answers  on 
both  computers. 

• Assure  concordance  of  all  RAYMODE 
versions  with  regard  to  the  physics 
used  • 

• Bring  RAYMODE  under  configuration 
management,  keeping  track  of  and  docu- 
menting all  upgrades,  cataloging  RAYMODE 
distribution,  and  providing  a feedback 
mechanism  for  identification  of  problems 
encountered,  errors  and  missing  but 
needed  features.  This  also  implies  a 
plan  for  upgrading  RAYMODE  or  otherwise 
altering  the  program.  It  is  particularly 
important  to  document  the  aspects  of 
RAYMODE  versions  which  differ  from  the 
basic  model  due  to  constraints  such  as 
achieving  a given  run  time  or  core  stor- 
age requirement.  It  is  also  required 
that  any  change  to  RAYMODE  undergo  test 
and  evaluation  before  being  distributed. 

(U)  No  model  evaluation  can  claim  to  be 
complete  and  such  is  the  case  here. 
There  are  some  particular  omissions  of 
test  case  scenarios  which  require  iden- 
tification: 

• Propagation  in  an  environment  char- 
acterized by  a sound  speed  profile  with 
a double  deep  sound  channel  such  as 
found  in  the  eastern  A.tlantic  Ocean, 

• No  analysis  was  performed  for  fre- 
quencies above  5 kilohertz  due  to  a lack 
of  experimental  propagation  loss  data 
with  supporting  environmental  data, 
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• No  analysis  for  shallow  water  scen- 
arios, once  again  due  to  a lack  of  ex- 
perimental data  (and  possibly  lack  of 
knowledge  of  the  physics  of  the  boundary 
interactions  where  their  effect  over- 
whelms the  effect  of  propagation  within 
the  water  medium), 

e Under-ice  propagation  was  not  exam- 
ined, and 

e The  deterioration  of  the  range  inde- 
pendent RAYMODE  X as  the  environment 
becomes  increasingly  range  dependent  was 
not  addressed  due  to  its  complexity. 

(U)  To  properly  evaluate  a model  such  as 
RAYMODE  X it  is  necessary  to  catalogue 
problem  runs  (i.e.,  those  that  produce 
clearly  invalid  answers),  the  frequency 
with  which  they  occur,  and  the  environ- 
mental and  acoustic  conditions  which 
accompany  their  occurrence.  Information 
of  this  type  is  most  valuable  for  prob- 
lem Identification  and  diagnosis  and  can 
mainly  be  obtained  from  user  feedback. 
The  importance  of  user  feedback  and  its 
solicitation  by  those  responsible  for 
configuration  management  of  the  RAYMODE 
model  cannot  be  overemphasized. 

(U)  Future  acoustic  experiments  should 
be  performed  with  model  evaluation  sup- 
port as  one  of  the  primary  objectives; 
this  has  implications  for  frequency  cov- 
erage, source/receiver  geometries,  data 
density  and  supporting  environmental 
measurements.  Attention  to  modes  of 
propagation  (e.g.,  surface  duct,  bottom 
bounce,  convergence  zone)  is  most  impor- 
tant to  model  evaluation  to  fill  many 
scenario  holidays. 
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1020  Balch  Boulevard 
Stennis  Space  Center  MS  39529-5005 

Subj : DECLASSIFICATION  OF  PARKA  I AND  PARKA  II  REPORTS 

Ref:  (a)  CNMOC  ltr  3140  Ser  5/1 1 0 of  1 2 Aug  97 

Enel : ( 1 ) Listing  of  Known  Classified  PARKA  Reports 

1 . In  response  to  reference  (a),  the  Chief  of  Naval  Operations  (N874)  has  reviewed  a number  of 
Pacific  Acoustic  Research  Kaneohe- Alaska  (PARKA)  Experiment  documents  and  has 
determined  that  all  PARKA  I and  PARKA  II  reports  may  be  declassified  and  marked  as 
follows: 

Classification  changed  to  UNCLASSIFIED  by  authority  of  Chief  of  Naval  Research 
letter  Ser  93/160,  10  Mar  99. 

DISTRIBUTION  STATEMENT  A:  Approved  for  public  release.  Distribution  is 

unlimited. 

2.  Enclosure  (1)  is  a listing  of  known  classified  PARKA  reports.  The  marking  on  those 
documents  should  be  changed  as  noted  in  paragraph  1 above.  When  other  PARKA  I and 
PARKA  II  reports  are  identified,  their  markings  should  be  changed  and  a copy  of  the  title 
page  and  a notation  of  how  many  pages  the  document  contained  should  be  provided  to  Chief 
of  Naval  Research  (ONR  93),  800  N.  Quincy  Street,  Arlington,  VA  22217-5660.  This  will 
enable  me  to  maintain  a master  list  of  downgraded  PARKA  reports. 

3.  Questions  may  be  directed  to  the  undersigned  on  (703)  696-4619,  DSN  426-4619. 


PEGGY  LAMBERT 
By  direction 

Copy  to: 

NUWC  Newport  Technical  Library  (Code  5441) 

NRL  Washington  (Mary  Templeman,  Code  5227) 

NRL  SSC  (Roger  Swanton,  Code  7031) 
l/DTIC  (Bill  Bush,  DTIC-OCQ) 


The  Acoustic  Model  Evaluation  Committee  (AMEC)  Reports,  Volume  2,  The  Evaluation  of  the  Fact 
PL9D  Transmission  Loss  Model,  Book  3,  Appendices  E-H,  September  1982,  NORDA-35-VOL-2-BK- 
1 3, 428  pages 
V (DTIC  # C034  020) 

The  Acoustic  Model  Evaluation  Committee  (AMEC)  Reports,  Volume  3,  The  Evaluation  of  the 
RAYMODE  X Propagation  Loss  Model,  Book  1,  September  1982,  NORDA-36-VOL-3-BK-1,  127 
pages 

(DTIC  # C034  021) 

The  Acoustic  Model  Evaluation  Committee  (AMEC)  Reports,  Volume  3,  TheEvaluation  of  the 
RAYMODE  X Propagation  Loss  Model,  Book  2,  Appendices  A-D,  September  1 982,  NORDA-36- 
VOL-3-BK-2,  324  pages 
(DTIC  # C034  022) 

The  Acoustic  Model  Evaluation  Committee  (AMEC)  Reports,  Volume  3,  The  Evaluation  of  the 
RAYMODE  X Propagation  Loss  Model,  Book  3,  Appendices  E-H,  September  1982,  NORDA-36- 
VOL-3-BK-3,  388  pages 
(DTIC  # C034  023)' 
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